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ABSTRACT 


The computer program COORDC generates a body-fitted curvi- 
linear coordinate system for corner geometry with or without 
corner fillets. It is assumed that at any given £, x remains 
constant; consequently the only variation is in y and z. It 
is also assumed that for all £'s in the physical plane the coordi- 
nate system in y-z plane is similar. This enables solution of 
coordinate system for one particular £ = 1 (x for £ = 1 is 
arbitrarily chosen to be 0.0) and the solution for all other 
£ plane can be easily specified once the coordinates in the 
physical plane on the line. l£££lMAX, ri = 1 , £ = 1 are specified. 
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List of Symbols 


IMAX 

Maximum number of points 

in 

K 

- direction 

JMAX 

Maximum number of points 

in 

n 

- direction 

KMAX 

Maximum number of points 

in 

c 

- direction 

J 

Jacobian 




P,Q,R 

Inhomogeneous terms used 

for 

coordinate attraction 

S,0,R 

Successive-Over-Relaxation iteration 

X,Y,Z 

Physical Coordinates 




a,3,Y 

Metric Coefficients 




K t n t C 

Transformed Coordinates 




r i' r 2 ,r 3 ,r 4 

Boundary contours in the physical plane 

r i*'V'V'V 

Boundary contours in the 

transformed plane 

£0 

Acceleration parameter for Gauss- 

-Seidel iteration 


Subscripts 

E> r T\ t £ r y , Z 

£C,nn,xx,yy etc. 
£n/xy/xz etc. 

Superscripts 

(s) 


Denotes first partial differentiation 
Denotes second partial differentiation 
Denotes cross partial differentiation 

Denotes the iteration number 


ii 



I. Numerical Generation of Body-Fitted 
Curvilinear Coordinate Systems 

The method for numerically generating boundary-fitted 
curvilinear coordinate systems is presented as applicable in 
a specialized 3-D corner with or without a fillet. The coor- 
dinate system is generated such that for a specified 5» x 
remains constant and for all 5's the solution in y-z plane 
stays similar. This is particularly beneficial because for 
5=1 plane (x is chosen to be 0.0 for 5=1) y-z solution 
can be obtained numerically by the 2-D formulation instead of 
3-D formulation of the body-fitted curvilinear coordinate system. 
The coordinates for all other 5 planes can be easily obtained 
once the physical plane coordinates on the line 5 = 1 » IMAX, 
n = 1, C = 1 are specified. In the present problem it is assumed 
that for a given 5 plane, right boundary and upper boundary have 
z and y constant respectively. In addition, on the left wall 
and the bottom wall, z and y respectively equal zero some dis- 
tance from the corner when fillet is present? in case of no fillet 
z and y on left and bottom wall respectively is equal to zero. 

The clustering of ti and £ lines is specified by the point dis- 
tribution on i) = 1 and z; = 1 line. 

Section A provides a brief discussion of the mathematical 
formulation and section B describes the numerical technique used 
to generate the coordinate system. Section C explains how the 
required initial guess is specified. Section D deals with handling 
of the x-coordinate and the metric coefficients. 




A. Mathematics of Transformation 


Consider the transformation of a simply connected region R, 

into a rectangular region as shown in Figure 1. We require that 

^ ^ ^ ^ 

^1' ^ 2 ' ^ 3 ' ^4 raa P ^ nto ^l • ^2 ' ^3 ' ^4 respectively. For 

identification purposes, region D is referred to as the physical 
* 

plane, and D as the transformed plane. In the physical plane 
it is assumed that x is constant for a specified £, z is constant 
on r 3 and y is constant on r^. If for £ = 1, x = 0.0 is assumed, 
then the transformation from x,y,z to 5 / 0 / C is defined as follows. 


~K~ 


1(x) - 

n 

= 

n (y,z) 

_ ? _ 


S (YrZ)_ 


The above equation is solved only for n and £ since 5 is 
explicitly known. Assuming that an inverse transformation 
exists, then 


( 1 ) 


y 

z 


y (n,c)" 
z (n ,e) 


( 2 ) 


where the jacobian J is as follows: 


J = y n z c ' y t 2 n 


(3) 


The generating elliptic system is chosen to be the inhomo- 
geneous Laplace equation. 


n + n = Q(n,c) 
yy zz ’ s ' 


(4a) 


C yy + ^zz = R(n ' ?) 


(4b) 



with the boundary conditions as follows: 


on r i 

T1 


(y»z) 



S 


_ *i 

(5a) 

on r 2 

n 


” n 2 



S 

1 

c 2 (y f z) 

(5b) 

on Tj 

T1 


~n 3 (y,zf 



s 


_ 

(5c) 

and on T. 

4 

~ r\~ 


r n “1 

4 



C 


c 4 (y^z) 

(5d) 


where ri 2 , are specified constants and r|^» C 2 , 

ti^/ C 4 are specified functions. 

The inhomogeneous terms Q and R are selected to control 
the spacing of n = constant and z, = constant lines in the physi- 
cal plane; several forms for the inhomogeneous terms can be used. 
In the present case a special form of Q and R terms was used based 
on the point distribution on r) = n 2 and z, — line. (For con- 
venience n 2 = = 1 and = JMAX and Z, = KMAX is chosen thus 

ensuring a field size of JMAX x KMAX for a given £ plane) . 

In order that the transformed plane contain a uniform grid 
of spacing equal to unity, the dependent and independent variables 



must be interchanged in equations (4a) and (4b) . The resulting 
equations after transformation are as follows: 


aY nn ‘ 26 y nc + Y y « 


az - 23 z + y z 

nn nc c? 


-J [y n Q(HfC)+ R(n,C)] 


-J [z n Q (n # c) + Z ? R(n,0] 


(6a) 

(6b) 


where , 


2 2 
a = y c + z 5 


= V V + z z 

y n n c 

2 2 
= y„ + z n 
n n 


j = v z — v z 


The boundary conditions transform as follows: 
* 


on r 

y 


a^Ti,^) 


z 


a 2 (n,^) 

r * 
on r 2 

”y ~ 


b 1 ( T l 2 , ?) 


z 


vv° 

* 

on r 3 

”y " 


"=1 (,1 ' C 3 S 


z 


o 2 (n,c 3 ) 

* 

and on r. 

4 

”y " 


"^1 (^4 f ^) 


z 

L — 


d 2 (n 4 r?) 

The functions a^, 

a 2 ' b l / 


(7a) 


(7b) 


(7c) 


(7d) 


A 



by known contours r F^, r^. If the coordinate point 

locations on these contours is fixed then dirchlet type boundary 

conditions can be used. On the other hand if the coordinate 

lines on these contours is to be made parallel or normal to 

* 

some line then Neumann type conditions can be used. On 
and dirchlet type conditions were used. On and 
Neumann type conditions were used. Thus the following expressions 

-k k 

are obtained on T_ and T. : 


on r 3 


y(J,KMAX) = y ( J,KMAX-1) (8a) 

z(J,KMAX) = constant z^ (8b) 

1<J<JMAX 

* 

and on F. 

4 

z(JMAX,K) = z ( JMAX-1 ,K) (8c) 

y(JMAX,K) = constant y^ (8d) 


1< K<_KMAX 

The last aspect to be discussed is how Q and R are computed. 
Q is computed based on the y variation of points on and R is 
computed based on the z variation of points on r 2 « The following 
discussion will make this point clearer. 

If one assumes on that coordinate points are specified 
such that the y derivatives with respect to £ are zero then fro™ 


(6a) , 

Q(n,C) 



(9a) 


Also, Q is not a function of C. Thus, only one one-dimensional 
array of size JMAX is sufficient to store Q. Similarly for R we 
obtain. 


/ 


-r- \ 



Again R is not a function of r|. A one-dimensional array 
of size KMAX is sufficient to store R. 


B. Numerical Solution Technique 

Based on the above discussion, the following equations must 
be numerically solved: 

o Y nn - 2 0 y n; + y y 5? = - J 2 [y^ Q(n) + R(?) ] (10a) 

a z - 2 $ z + y 2 = - J 2 [z Q(n) + z R(?) ] (10b) 

no n5 ss n ? 

where 

Q(n) = - on (11a) 

y n 

and 

R(?) - - on T (lib) 

z 

C 

The first and second order derivaties in (10) and (11) are 
written in central difference form. The expressions for Q and R 
are as follows: 


y (J+l, 1) - 2*y(J,l) + y ( J-l , 1) 
(0 . 5* (y ( J+l , 1) - y(J-l,l) ) )**3 

2< J< JMAX-1 


(12a) 


and 


R(K) 


z (1,K+1) - 2*z (1 ,K) + z (1 ,K-1) 
( 0 . 5* (z ( 1 ,K+1) - z (1 ,K-1) ) ) **3 

2< K< KMAX-1 


(12b) 


Equations (10a) and (10b) are solved by successive over 



relaxation (S.O.R) technique. The intemediate value for y and 
z are as follows: 


Y(J,K) = [a (y (J+1,K) + y (J-1,K) ) - 2 8 

+ y(y (J/K+l) + y (J ,K-1) ) + J 2 (y Q(J) 

+ y ? R(K)) ] / [ 2 (a+y) ] (13a) 

where 

y n? = 0.25 ty(J+l,K+l) - y (J+1,K-1) - y (J-l/K+1) 

+y(J-l,K-l)] (13b) 

y n =0.5 [y (J+1,K) - y '.J-1,K) ] (13c) 

y c =0.5 [y (J /K+l) - y ( J ,K-1) ] (13d) 

a = 0.25 [ (y (J,K+1) - y (J ,K-1) ) 2 (Z(.T f K+l) - z(J,K-l)) 2 ] 

8 = 0.25 [ (y (J+l, K) - y(J-l f K)) (y(J,K+l) - y(J,K-l)) 

+ (z(J+l,K) - z(J-l,K) ) (z(J,K+l) - z (J ,K-1) ) ] 
y = 0.25 [ (y (J+l ,K) - y (J-1,K) ) 2 + z(J+l,K) - z(J-l,K)) 2 ] 

J = 0.25 [ (y (J+1,K) - y(J-l,K)) (z(J,K+l) - z ( J ,K-1) ) 

- (y (J f K+l) - y(J,K-l)) (z(J+l,K) -z(J-l,K))] 


similarly an expression for z can be written out. For each grid 


(13e) 

(13f ) 
(13g) 

(13h) 



point two finite difference equations (one for y and one for z) 
are obtained. The subscripts J and K have a range as follows: 

J = 2,3/4, , JMAX-1 

K = 2,3,4, , KMAX-1 

This results in a set of (JMAX-2) (KMAX-2) . The above equa- 
tions are solved simultaneously using S.O.R. For a set of non- 
linear differential equations for some variable f, S.O.R. itera- 
tion can be written as follows: 

f ( s +D _ w .f + (i _ W ) f (s) (14) 

where superscripts denote the iteration number, w is the 
acceleration parameter, and the latest values are used in the 
difference equations to solve f . In the present problem f 
denotes y or z. 

Since on a nd * 4 , Neumann type boundary conditions are 
applied, these conditions are implemented as follows: 
on r 3 

y ( J,KMAX) = y(J, KMAX-1) (13a) 

2<J<JMAX 

on r . 

4 

z (JMAX,K) = z (JMAX-1, K) (15b) 

2<K<KMAX 

Since on z = constant and on y = constant, equations 
(15a) and (15b) ensure that coordinate lines on and are 

orthogonal even though the coordinate system is not necessarily 
orthogonal everywhere. 


R 



C. Initial Guess and Boundary Point Specification 

The input to the program consists of points specification 
on r and boundary. boundary is such that the variation 

of z close to the JMAX point is zero. Similarly on ^ the y 
variation close to KMAX point is zero. Points on and ^ are 
specified implied that (y(J,l), z(J,l)} 1<J<JMAX and {y(l,K), 
z(l,K)} 1<_K<_KMAX are known. Then for the initial guess we have/ 

z(J,KMAX) = z (1,KMAX) , 2<J<JMAX (16a) 

y (JMAX/K) = y (JMAX, 1) , 2<K<KMAX (16b) 

also, 

y(J,K) = y (1,K) + (y(JMAX,K) - y(l,K)) * ((y(J,l) 

-y (1 , 1) ) / (y ( JMAX, 1) - y (1,1) ) ) (16c) 

z (J,K) = z(J,l) + (z(J,KMAX) - z(J,l) ) * ( (z(l,K) 

-z(l,l))/(z (1,KMAX) - z(l,l))) 

2<_J<_JMAX-1 and 2<K<KMAX-1 

and 

y (J, KMAX) = y (J,KMAX-1) , < J<JMAX-1 

Z (JMAX , K) = z ( JMAX-1 ,K) , <_K<_KMAX-1 

D. x-Coordinate and Metric Coefficients 

As discussed before, the y and z variation is computed 
only for one plane 5=1; also x is constant for any specified 
5-plane. Thus if x,y,z variation on 1<_5^_IMAX, n = 1 and t, = 1 
is specified then coordinates of any point can be easily speci- 
fied. For example say for 5 = II, 


(16d) 


(16e) 

(16f) 



x(Il,l,l) = xl (17a) 
y(n,U) = yi (17b) 
2(11,1,1) = zl (17c) 

Then for any point II, J, K in £ = II plane we have, 

x (II, J,K) = xl (17d) 
y(Il,J,K) = y (11,1,1) - y (1,1,1) + y(l,J,K) (17e) 
z (II, J ,K) = z(Il,l,l) - z(l,l,l) + z (1, J,K) (17f ) 


If the coordinate system on 5 = 1 plane is continuous 
in y and z, and if on 5 = 1, IMAX, n = 1, ? = 1 line the var- 
iation in x, y, z is continuous then the solution given by 
above equations will be continuous. This is ensured by having 
continuous y and z variation on and and on line 5=1/ 
IMAX, n = 1, 5=1. 

A brief discussion on metric coefficients can now be 
included. The nine metric coefficients 3. . are given by: 


11 

<y n 

z c 

- y 5 

z n )/ J 

12 

(y 5 

z c 

- y s 

Z ? )/J 

13 

(y c 

Z 

n 

" y n 

z c )/J 

21 = 

(x c 

z 

n 

- X 

n 

Z^)/J 

22 

(x f 

z ? 


Z ? )/J 


(13a) 

(13b) 

(13c) 

(lOd) 


(18e) 



Since x is constant for any specified £ plane, we have 
x^ = = 0 which considerably simplifies the above equations. 

Also, we can ensure that y^ = 0, z^ = 0 if v and z are specified 
to be constant on 1<_£<IMAX, n = 1, 5 = 1 line. In such a case 
the metric coefficients become. 


*11 


(y n 2 ? ' y 5 V /J 

(19a) 

*12 

= 

0.0 

(19b) 

*13 

= 

0.0 

(19c) 

*21 

= 

0.0 

(13d) 

*22 

= 

X £ V J 

(19e) 

*23 

= 

- X- z /J 

£ rr 

(l?f ) 

*31 

= 

0.0 

(ISg) 

*32 

= 

" X £ y 5 /J 

(19h) 

*33 


X E V J 

(19i) 

J 

= 

V y n 2 ? - '(V 

(19? ) 


E. Concluding Remarks 

A method for generating 3-D body-fitted coordinate systen 
for corner flows has been presented. The computer code to 
generate coordinate systems using equations 12, 13, 14, 15, 1C 
is presented in Chapter 3. The coverge'’ ■'•'^ordinate system is 
stored on a disk file which becomes a part of the input to the 
Navier Stokes code. Successful coordinate 


svstems for coners 



with or without fillet and for a 3-D flat plate have been 
generated (the program can generate suitable coordinate 
systems for a 3-D flat plate problem) . The compressible 
laminar Navier Stokes solution for the above mentioned 
problems at low Reynolds number have already been obtained. 
Effort is now underway to obtain solution for higher Reynolds 
number with or without turbulence. 







II. Computer Program 
Listing of the Computer Program 
The program listing follows this page 



REATICN RUN 


CAROS ENCOUNTERED IN INPUT 


PAGE X 




UPDATE 1*2-452* 


tOECK CORD1 

PROGRAM COORDC ( INPUT* OUTPUT , SOL NX X, TAPE l«$OLNXX, TAPES* IN PUT* 

1 TAPE6*0UTPUT) 

C + + 

£*****+****+*********************+**+**♦+++**♦+*%+*+*+♦+*+++♦+♦**+**♦ 


C + + 
C + + 
C+ + 
C + + 
C + + 
C + + 
C + + 


* 

CORNER COORDINATE GENERATION PROGRAM BASED ON METHOD ♦ 

DEVELOPED BY THOMPSON, THAMES > AND MASTIN OF MISSISSIPPI * 
STATE UNIVERSITY* * 

♦ 

FOR FURTHER ENQUIRIES CONTACT l ♦ 

♦ 


C + + 
C + + 
C + + 
C + + 

C + + 
C + + 
C + + 


DR* JULIUS HARRIS 
DR. DILIP KUMAR 

NASA - LANGLEY RESEARCH CENTER 
MAIL STOP 163 
HAMPTON, V A* 23665 
PHONE* 804-827-3696 
FTS CODE 928-1110 


+ 

+ 

♦ 

♦ 

♦ 

♦ 

* 


C + + 

C + + 

£***+♦*** ***** 


♦ 

tt*+++*t*+*+*++++**+++**+++++++4++ **+++++** ♦♦+♦♦♦+*♦♦♦♦ 


C + + 

C + + 

C + + 

C + + 

C + + 
C+ + 

C + + 

C + + 

C + + 

C + + 

C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 

c** 


THIS PROGRAM GENERATES COORDINATE SYSTEM FOR A 2-D CORNER 
WITH OR WITHOUT FILLER. THE METHOD USED IS BASED ON METHOD 
DEVELOPED BY THOMPSON/ THAMES/ MASTIN AND OTHERS AT MISS. 
STATE UNIVERSITY. THE CORNER IS SUCH THAT EXCEPT IN THE 
FILLER (WHEN PRESENT) L.H.S, BOUNDARY HAS 2 • 0 AND BOTTOM 
SIDE BOUNDARY HAS Y - 0. ALSO R.H.S. BOUNOARY ALWAYS HAS 
2 ■ CONSTANT AND UPPER BOUNDARY HAS Y « CONSTANT. THE 
PROGRAM GENERATES Q AND R SUCH THAT COORDINATE CLUSTERING 
IS PROPORTIONAL TO COORDINATE POINTS DISTRIBUTION ON L.H.S. 
AND BOTTOM BOUNDARY, ON R.H.S. AND TOP BOUNDARIES NEUMANN 
TYPE BOUNDARY CONDITIONS ARE USED. THIS ENSURES THAT 
COORDINATE LINES ON R.H.S. AND UPPER BOUNDARY ARE PARALLEL 
TO SIDES. 


CARD INPUT 

CARDll SPECIFIES FIELD SIZE AND PRINT/ PLOT FLAGS 

JHAX,KHAX,IPRU>IPRT2f IPLTI/IPLT2/ITERMX/R1/R2/R3 
( FORMAT 7 15 » 5X » 3F10.0 ) 


79/06/26* 

13*30*37 

CORDl 

1 

com 

2 

com 

3 

com 

4 

com 

5 

com 

6 

com 

7 

com 

8 

com 

9 

com 

10 

com 

11 

com 

12 

com 

13 

com 

14 

com 

15 

com 

16 

com 

17 

CORDl 

18 

C0RD1 

. 19 

CORDl 

20 

C0R01 

21 

CORDl 

22 

CORDl 

23 

CORDl 

24 

CORDl 

25 

CORDl 

26 

CORDl 

27 

CORDl 

20 

CORDl 

29 

CORDl 

30 

CORDl 

31 

CORDl 

32 

CORDl 

33 

CORDl 

34 

CORDl 

35 

CORDl 

36 

CORDl 

37 ' 

CORDl 

30 

CORDl 

39 

CORDl 

40 

CORDl 

41 

CORDl 

42 

CORDl 

43 

CORDl 

44 




: AT ION RUN 

C++ 
C + + 
C + + 
C++ 
C + + 
C++ 
C+ + 
C + + 
C + + 
C++ 

c+* 

C + + 
C + + 
c»* 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C+ + 
C+ + 
C + + 
C + + 
C + + 
C+ + 
C + + 
C + + 
C+ + 
C + + 
C+ + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 
C + + 


CAROS ENCOUNTERED IN INPUT 


UPDATE 1.2-452. 


JMAX I 
KMAX t 
IPRTI i 


IPRT2 i 
IPLT1 x 
IPLT2 : 
ITERMX t 
R1 I 

R2 l 

R3 l 


NO OF NODES IN Y-DIRECTION 
NO OF NODES IN Z-DIRECTION 
INITIAL GUESS IS TO BE PRINTED OR NOT 
IF IPRTI • 1 THEN PARTIALLY CONVERGED 
SOLUTION IS PRINTED EVERY 50 ITER# 
FINAL SOLUTION TO BE PRINTED OR NOT 
INITIAL GUESS TO BE PLOTTED ON NOT 
FINAL SOLUTION TO BE PLOTTED OR NOT 
MAXIMUM NO. OF ITERATIONS ALLOWED 
GAUSS SEIDELL ITERATION PARAMETER 
(USUALLY 0*6) 

CONVERGENCE CRITERION IN Y-OIRECTION ' 
(USUALLY 0.000005 TO 0.00001) 
CONVERGENCE CRITERION IN Z-DIRECTION 
(USUALLY SAME AS R2) 


NOTE x FOR ALL FLAGS 0 MEANS NO# 1 MEANS YES 


CARD2I CARDS TO SPECIFY POINT DISTRIBUTION ON L.H.S. 
BOUNDARY (Y-OIRECTION) NO OF CARDS * JMAX 
Y/Z (FORMAT 2F10.0) 

Y l Y COORDINATE 

Z x l COORDINATE 

CARD3: CARDS TO SPECIFY POINT DISTRIBUTION ON BOTTOM 

BOUNDARY (Z-DIRECTION). NO. OF CARDS ■ KMAX. 

YiZ (FORMAT 2F10.0) SEE CARO TYPE 2) 

-CARO 4* SPECIFIES NO. OF NODES IN X-DIRECTION. 

IMAX ( FORMAT ( 15 ) 

IMAX I NO. OF POINTS IN X-DIRECTION. 

CARDS * CARDS TO SPECIFY NO. OF POINTS IN X - OIRECTION. 

X> Yl/ Zl (FORMAT 3F10.0) 

X t X-COORDINATE ON J-1,K»1 LINE ( X-0 IRECTION) 

Y1 l Y-COORDINATE ON J-1*K»1 LINE ( X-DIRECTION) 

ZI I Z-COOROINATE ON J-l/K-l LINE (X-DIRECTION) 

CARD6* PRXNT/PLOT HE AOINGS • 

HED1(4) (FORMAT <i A10 ) 

HED1 I 40 BYTES OF PLOT HEADING 
CARp7l ADDITIONAL PRINT PLOT HEADINGS. 

HEDZm (FORMAT 4A1Q) 



79/06/26. 13 

.30. 37 < 

C0R01 

45 

C0R01 

46 

C0R01 

47 

C0RD1 

48 

C0RD1 

49 

C0RD1 

50 

C0RD1 

51 

C0RD1 

52 

C0RD1 

53 

C0RD1 

54 

C0RD1 

55 

C0RD1 

56 

CORD! 

57 

C0RD1 

58 

C0RD1 

59 

C0RD1 

60 - 

C0RD1 

61 

C0RD1 

62 

C0RD1 

63 

C0RD1 

64 

C0RD1 

65 

C0RD1 

66 

C0RD1 

67 

C0RD1 

68 

C0RD1 

69 

C0RD1 

70 

C0RD1 

71 

C0RD1 

72 

C0R01 

73 

CDRD1 • 

74 

C0RD1 

75 

cordi 

76 

C0RD1 

77 

CORDI 

78 

CORDI 

79 

CORDI 

80 , 

CORDI 

81 

CORDI 

82 

CORDI 

83 

CORDI 

84 

CORDI 

05 

CORDI 

86 

CORDI . 

87 

CORDI 

88 




i • 


RUN 

CAROS ENCOUNTERED IN INPUT 

UPOATE 1.2-452. 

79/06/26, 

13,30,37# ; 

C++ 

HE02 1 40 BYTES OF ADDITIONAL PLOT HEADING. 

C0RD1 

89 

C++ 


• 

CORD! 

90 

C++ 

SAMPLE RUNSTREAM. 


C0RD1 

91 

C + + 



C0RD1 

92 

C+ + 



C0RD1 

93 

C + + 

KVIS2# T777# CM70000, BIN 34 

D KUHAR 

C0RD1 

94 

C + + 

U$ER#479Q19C* 


C0RD1 

95 

C + + 

CHARGE#100718#LRC, 


C0RD1 

96 

C + + 

GET ( 0LDPL»KMRCD1/UN*375732N) • 


com 

97 

C + + 

UPDATE(F#C«VISFL1), 


C0RD1 

98 

C + + 

C0PYBR#INPUT#EDTFT1, 

. 

C0RD1 

99 

C+ + 

REWIND# VI SF L 1# EDTFT1, 


C0RD1 

100 

C + + 

EDIT#VISFL1#I-E0TFT1. 


C0RD1 

101 

C + + 

REWIND# VISFL1 , 


COROl 

102 

C + + 

FTN(0PT-1#I*VISFL1), 


CORD! 

103 

C + + 

ATTACH#LRCGO$F/UN-LIBRARY#NA, 


C0RD1 

104 

C + + 

LDSET#LlB-LRCGO$F#PRE$EWERO, 


CORDI 

105 

C + + 

LGO. 


CQRDl 

106 

C + + 

PLOT • VARI AN 


C0RD1 

107 

C + + 

REWIND# 50LNXX « 


C0RD1 

108 

C+ + 

C0PY8R#S0LNXX#FILETC, 


C0RD1 

109 

C+ + 

REWIND#FILETC. 


CORDI 

110 

C + + 

SAVE#FILETCt 


C0RD1 

111 

C+ + 

‘7/8/9 


CORDI 

112 

C + + 

♦IDENT 0UM1 


CORDI 

113 

C + + 

♦INSERT COROl, 3 


CORDI 

114 

C + + 

C + + 


CORDI 

115 

C + + 

7/8/9 


CORDI 

116 

C + + 

R $1/ HXLL /#/ 30/) + 


CORDI 

117 

C + + 

RSI/ MXLL /# / 23/ J * 


COROl 

118 

C + + 

RSi/ MXLL/#/ 23/J+ 


CORDI 

119 

C++ • 

F.NO 


CORDI 

120 

C + + 

7/B/9 


CORDI • 

121 

C + + 

INPUT DATA, SEE EXPLANATION ABOVE# 


CORDI 

122 

C + + 

6/7/B/9 


COROl 

123 , 

C + + 



COROl 

124 

C + * 



COROl 

125 

C + + 

. 


CORDI 

126 

c+* 

REWIND# VISFLI, 


CORDI 

127 

C + + 

NOTEl AFTER FIRST 7/8/9 CARD 3 CARDS ARE INPUT TO 


CORDI 

128 

C + + 

SATISFY INPUT FOR UPDATE UTILITY, 


CORDI 

129 

C + + 



. CORDI 

130 

C + + 

AFTER 2ND 7/8/9 CARD EDIT CARDS ARE INPUT, 

THE 

CORDI 

131 

CM 

PROGRAM IS WRITTEN SO THAT MAXIMUM NUMBER OF 


CORDI 

132 


:reation~k un Tttorwmmzms in inpui — — prorrr u.a-re zT 


c** 
c** 
c** 
c ♦* 

C ** 

c** 

C** 

c*+ 
c + * 
c** 
c** 
c*+ 
c** 
c** 
c** 
c** 
c** 
c** 
c** 
c** 
c** 
c** 


POINTS IN X,Y,Z DIRECTION CAN 8E CHANCED BY SIMPLE 
INPUT. IN THE ABOVE CASE FIELD SIZE OF 30/23,23 
IS ASSUMED. IN THE THREE CARDS AFTER RS«/ AND MXU 
I/J/K SHOULD BE INPUT RESPECTIVELY. THIS COULD NOT 
BE SHOWN ABOVE BECAUSE TEXT EDITOR CHANGES TEXT IN 
COMHENTS CARDS ALSO. 

IT IS ASSUMED THAT PROGRAM RESIOES ON FILE KMRCD1 
UNDER USER NO. OF 375732N AND USER 479019C HAS 
PERMISSION TO USE IT. FURTHERMORE, IT IS ALSO ASSUMED 
THAT THE COORDINATE SYSTEM SOLUTION IS TO BE STORED 
ON FILE FILETC. THIS NAME CAN BE CHANGED BY 
CHANGING THREE CARDS TOWAROS THE END OF THE RUNSTREAH, 

THE PLOTS ON CALCOMP CAN BE 08TAINE0 BY REPLACING 
PLOT.VARIAH CARO BY FOLLOWING CARDS. 
PL0T.CALP0ST,11(X0»1.0, Y0-0.2) 

CONT. //PAPER 00, 

CONT. LEROY PEN SIZE 0.3 

CONT. MULTIPLE PLOT MOOE 

CONT. THESE PLOTS ARE FOR REPORTS// 


C** 


c** 


c** 

c** 

c** 


1 

1 


COMMON /BXXX / Y ( JMXLL/KMXLL ) / Z ( JMXLL, KMXLU/Q( JMXLL), RCKHXLL), 
X(IMXLL)/Y1(IMXLL)/Z1(IMXLL) 

COMMON /AXXX/ ILIMIT, JLIMIT, KLIMIT, IKAX, JHAX,KMAX,IPRT1,IPRT2, 
IPLT1>IPLT2/ITERHX»R1,R2/R3/I1/HE01(A)/HED2(A) 


DATA KKK, YERRMX, ZERRMX,I1,IYL0CJ, IYLOCK, IZLOCJ/IZLOCK /O, 2*0., 5*0/ 
ILIMIT • IMXLL 
JLIMIT • JHXl.L 
KLIMIT * KMXLL 


INITIALIZE 0, R, Y AND Z TO ZERO, 


DO 10 J » 1, JLIKIT 
ON) • 0,0 
10 CONTINUE 

00 20 K • 1/KUHIT 

R (K ) « 0.0 
20 CONTINUE 

00 30 J • 1, JLIMIT 
00 30 K • 1,KLIMIT 
Y(J,K) « 0.0 


7970 677 6VT373 OT 37T~”' w ~PAG E 


C0RD1 

133 


C0RD1 

13* 


C0RD1 

135 


C0RD1 

136 


CORD! 

137 


C0RD1 

138 


CflRDl 

139 


COROl 

1*0 


C0RD1 

1*1 


COROl 

1*2 


C0RD1 

1*3 


COROl 

1** 


COROl 

1*5 


COROl 

1*6 


COROl 

1*7 


COROl 

1*8 


COROl 

1*9 


COROl 

150 


COROl 

151 


COROl 

152 


COROl 

153 


COROl 

15* 


COROl 

155 


COROl 

156 


COROl 

157 


C0RD1 

158 


COROl 

159 


C0RD1 

160 


C0RD1 

161 


C0RD1 

162 


COROl 

163 


COROl 

16* 


COROl 

165 


COROl 

166 


COROl 

167 


COROl 

168 


COROl 

169 


COROl 

170 • 


CORD! 

171 


COROl 

172 


COROl 

173 


COROl 

17* 


COROl 

175 


COROl 

176 



CREATION RUN 


CAR0S E NCOUNTERED IN INPUT UP QATE 1.2-452, 79 /06/2 6. 13 .30.37. 


.PACE 5 


3 



CAROS ENCOUNTERED IN INPUT 


UPDATE 1.2-452. 


79/06/26. 13.30.37 


CREATION RUN 



Z(J,K) • 0.0 


CORDl 

177 

30 

CONTINUE 


CORDl 

178 

C** 



CQR01 

179 

c ** 

READ INPUT 


C0R01 

180 

C ** 



CORDl 

181 


READt 5,1000) JMAX»KMAX,IPRT1,I PRT2, IPLT1, IPLT2, ITERHX.R1.R2.R3 

com 

182 


IFtJMAX.GT.JLIMIT.OR.KMAX.GT.KLIMIT) GO TO 50 


CORDl 

183 


GO TO 60 


CORDl 

184 

C** 



CORDl 

185 

50 

WRITE (6. 1020) JMAX.KMAX.JUMIT.KIIHIT 


CORDl 

186 


CALL EXIT 


CORDl 

187 

C + * 



CORDl 

188 

60 

CONTINUE 


CORDl 

189 


RE AD (5. 1030) C C Y( J.D.IU.Dl.M. JHAX) 


CORDl 

190 


READ( 5,1030) C ( Y (l.K ), Z ( l.K ) ), K-l.KHAX ) 


CORDl 

191 

c** 



CORDl 

192 

c** 

READ IN POINTS ALONG J « 1, K • 1 LINE (X-DIRECTION) 

CORDl 

193 

c* + 



' CORDl 

194 


READ! 5,1000 ) IMAX 


CORDl 

195 


IFt IHAX.LE.ILIMIT) GO TO 65 


CORDl 

196 


WRITE l 6, 1090 ) IMAX.ILIHIT 


CORDl 

197 


STOP 


CORDl 

198 

65 

CONTINUE 


CORDl 

199 




CORDl 

200 


READ ( 5,1100) (<Xm,YlCI),Zl(I)),I»l,IHAX) 


CORDl 

201 

C + * 



CORDl 

202 

C** 

SPECIFY INITIAL GUESS ON THE WHOLE FIELD 


CORDl 

203 

c+* 



CORDl 

204 


DO 70 J • 2, JHAX 


CORDl 

205 


ZIJ.KMAX) « ZCl.KHAX) 


CORDl 

206 

70 

CONTINUE 


CORDl 

207 


DO 80 K » 2.KHAX 


CORDl 

208 


Y(JHAX.K) * Y t JHAX, 1 1 


CORDl 

209 

‘80 

CONTINUE 


CORDl 

210 

C** 



CORDl 

211 


J1 « JHAX - 1 


CORDl 

212 


K1 • KHAX - 1 


CORDl 

213 


DO 100 J • 2, J1 


CORDl 

214 


DO 100 K * ?,K1 


CORDl 

213 


Y(J.K) ■ (Y( JHAX.K) - T(1,K)) ♦ (tY(J.l) - YCl/lU 

/ ' 

CORDl 

216 


1 (YUMAX.l) - Ytl.l ) )) ♦ Y(l.K) 


CORDl 

217 


Z(J,K) ■ (Z(J.KHAX) - Z(J.l)) ♦ l(Itl,KI - Z ( 1# 1 ) ) 

/ 

CORDl 

218 


1 (Z(l.KHAX) - Z(l,l)M ♦ Z(J,1) 


CORDl 

219 

100 

CONTINUE 


CORDl 

no 



CREATION RUN 


CAROS ENCOUNTERED IN INPUT UPDATE 1.2-A52. 79/06/26. 13.30.37. 


DO 102 J • 2,J1 

Y { J.KMAX ) • Y(J.KHAX-l) 

102 CONTINUE 

DO 10A K » 2>K1 

Z ( JMAX.K) - Z(JMAX-1,K> 

10A CONTINUE 
C** 

READ PRINTjPLOT HEADINGS. TWO CARDS ARE READ IN, 
c ** CALL PLOT INITIALIZING ROUTINE PSEUDO IF PLOT OPTIONS 

C** SET. 

C** 


READ(5,1080> (HED1( I ) > 1*1, A J 
READ(5>1080) (HE02< I ) t I«1 , A1 
IF(IPLT1,E0,0,AN0.IPLT2.E0.0) C-0 TO 110 
CALL PSEUDO 
CALL FONTS (1 ) 

110 CONTINUE 
C** 

c** SEE IF INITIAL GUESS NEEOS TO BE PRINTED 

C** 

IF{ IPRT1.E0.1) CALL I PRTC (l/KKK; YERRMX/ZERRHX* IYLOCJ#IYLOCK, 

1 I ZLOC J » I ZLOCK ) 

C* + 

C*» SEE IF INITIAL GUESS NEED TO BE PLOTTED 

C** 

IF(IPLTI.EO.I) CALL IPLTC (1;KKK; YERRMX/ZERRMX; IYLOC J» IYLOCK» 
l IZLOCJ.IZLOCK) 

C** 


C + + 
C + + 
C+ + 
C++ 


.COMPUTE 0 AND R. FIRST YETA OR ZTAU 1$ COMPUTED# 
Q ■- YETAETA / < YETA+ YETA+YETA ) 

R *- ZTAUTAU / ( ZTAU + ZT AU + ZTAU ) 


DO 1Z0 0 
A1 » 


0(J) 

120 CONTINUE 
C + + 


- 2/J1 

0.5 * (YCJ+1,1) - YU-hl)) 

■ - (Y ( J + lil ) - 2.0+Y( J# 1) ♦ YCJ-1#!)) / ( Al+Al+Al) 


DO 140 K ■ 2,K1 

u • 0.5 * < zci/K+i j - zu,K-m 

R(K) • - (ZC1/K+1) - 2*0+Z(l>K) ♦ Z(lfK-ll) / (AltAl+Al) 
140 CONTINUE 
C + + 

'C++ START COMPUTING SOLUTION. 

C+ + 


C0R01 

221 

CORDX - 

222 

CQRDi 

223 

C0RD1 

224 

C0RD1 

225 

C0RD1 

226 

CQRDI 

227 

C0RD1 

228 

C0RD1 

229 

C0RD1 

230 

C0RD1 

231 

C0RD1 

232 

C0RD1 

233 

C0RD1 

234 

C0RD1 

235 

C0RD1 

236 

C0RD1 

237 

C0RD1 

238 

CQRDI 

239 

C0RD1 

240 

C0RD1 

241 

C0RD1 

242 

C0RD1 

243 

CQRDI 

244 

CORDl 

245 

CQRDI 

246 

CORDl 

247 

CORDl 

248 

CORDl 

249 

CORDl 

250 

CORDl 

251 

CORDl 

252 

CORDl 

253 

CORDl 

254 

CORDl 

255 

CORDl 

256 • 

CORDl 

257 

CORDl 

258 

CORDl 

259 

CORDl 

260 

CORDl 

261 

CORDl 

262 

CORDl 

263 

com 

264 # 


PAGE 


CREATION RUN 

CAROS ENCOUNTERED IN INPUT UPDATE 1,2-452. 

79/06/26* 

13*30*37 

KKK 

* 0 

C0RD1 

265 

I ERR • 0 

CORDl 

266 

DO 

500 11 - If ITERMX 

C0RD1 

267 


YERRMX » 0.0 • 

C0R01 * 

260 


ZERRMX • 0.0 

C0RD1 

269 


CALL CALCOR(YERRHX/ZERRHXfIYLOCJfIYLOCKfIZLOCJ/IZLOCK) 

C0RD1 

270 

C + + 


C0RD1 

271 

C++ 

CHECK IF MAX ERROR TO BE PRINTED OR NOT 

C0RD1 

272 

c+* 


C0RD1 

273 


IERR « I ERR + 1 

C0RD1 

274 


IFIIERR.NE.50) GO TO 160 

C0RD1 

275 


IERR • 0 

C0RD1 

276 


WR I TE ( 6f 10A0 ) IIfIYLOCJfIYLOCKfYERRMXfIZLOCJfIZLOCKfZERRMX 

C0RD1 

277 

I Ft IPRT1. EO, 1 > CALL IPRTC (2 /KKKfYERRMXfZERRMXf IYLOC IYLOCKf 

C0RD1 

278 

1 

IZLOCJf IZLOCK) 

CORD! 

279 

160 

CONTINUE 

CORDl 

280 

C + + 


C0RD1 

201 

C + + 

CHECK FOR CONVERGENCE 

C0RD1 

282 

C + + 


CORDl 

283 


IF(YERRMX,GT«R2) GO TO -180 

CORDl 

284 


IFIZERRMX.GT.R3) GO TO 180 

CORDl 

285 

C + + 


CORDl 

286 

C + + 

SOLUTION CONVERGED, GET OUT OF THE LOOP. KKK • 1 

CORDl 

287 

C + + 

indicates solution converged 

CORDl 

288 

C + + 

•v. • 

C0PD1 

289 


KKK • 1 

CORDl 

290 


GO TO 520 

CORDl 

291 

C + + 


CORDl 

292 

180 

, CONTINUE 

CORDl 

293 

500 CONTINUE 

CORDl 

294 

520 CONTINUE 

CORDl 

295 

C + + 


CORDl 

296 

C + + 

see if converged or partially converged solution is to 

CORDl 

297 

C + + 

BE PLOTTED and PRINTED OR NOT 

CORDl 

290 

C + + 


CORDl 

299 

IF(IPRT2.E0,1) CALL IPRTC (2fKKKfYERRMXfZERRMX/ IYLOC J» IYLOCKf 

CORDl 

300 # 

1 

IZLOCJ/IZLOCK) 

CORDl 

301 

C + + 


CORDl 

302 

IFt IPLT2.EQ.1) CALL IPLTC C2>KKK>YERRMX> ZERRHX# IYLOC J/ IYLOCKf 

CORDl 

303 

1 

I ZLGC J ; I ZLOCK) 

CORDl 

304 

C + + 


CORDl 

305 

C + + 

STORE CONVERGED SOLUTION 

CORDl 

306 

C + + 


CORDl 

307 

IF(KKK.EOtO) GO TO 540 

CORDl 

308 



CREATION RUN 


CAROS ENCOUNTERED IN INPUT UPDATE 1*2-452* 


€*♦ 


WRITE (1/1110) ILIMIT/ JLIMIT/KLIHIT/ IHAX/ JHAX/KMAX/ 

1 <(HED1(I)/HED2(I))/I«1/A) 

WRITE (1/1120) ( ( (Y( J/K)/Z(J»K) )/K»l/KLIMIT)/J«l/JLIHIT) 
WRITE (1/1120) ((X(1)/Y1(I)/Z1(I))/I«1/ I LIMIT) 

GO TO 560 


SOLUTION 010 NOT CONVERGE. WRITE A MESSAGE 

C + * 

5 AO CONTINUE 

WRITE (6/ 1060) YERRMX/IYLOCJ/ IYIOCK/ZERRMX/ IZLOCJ/ IZL0CK/R2/R3/ 
1 ITERMX 

C** 

C** ALL FUNCTIONS COMPLETE 

C + * 


560 CONTINUE 
STOP 

C* + 


C** FORMAT STATEMENTS 

C+* 


1000 FORMAT(7I5/5X/3F10.0) 

1020 F0RHAT("1"////5X/"*** NO OF POINTS IN Y OR Z GREATER THAN"/ 

1 " MAXIMUM ALL0WED» n //5X/"*t* POINTS INPUT • "/2I5///5X/ 

2 " MAXIMUM ALLOWED • "/2I5/ //5X."*** RUN ABORTED") 

1030 FORMAT ( 2F10.0 ) 

10A0 FORMAT ("0"/5X/» ITER NO./ I YLOC J/ IYLOCK/ YERRMX/ IZLOCJ/ IZLOCK/"/ 

1 "ZERRMX - "/3I5/E15.6/2I5/E15.6 ) 

C** 

1060 FORMAT ( "l"////5X/"+** SOLUTION DID NOT CONVERGE. SOLUTION"/ 

1 . " NOT STORED ON UNIT 1"//5X/"*** MAX ERROR IN Y ■ "/E15.6/ 

2 " OCCUREO AT LOCATION • "/ 15 » "»"» I5//5X/ •'*♦* MAX ERROR IN " 

3 »"Z . "/ El 5 . 6/ " OCCUREO AT LOCATION • "/ I5/"/"/I5//5X/ 

A "*** MAX ERROR ALLOWED IN Y AND Z « "/ 2E15.6//5X/ "*♦* NO."/ 

5 " OF ITERATIONS PERFORMED ■ »»I5) 

1080 FORMAT ( AA10 ) 

1090 FORMAT ("1"/ ///5X/"*** NO. OF POINTS IN X-DIRECTION GREATER THAN"/ 

1 " MAXIMUM ALLOWEO ."// 5X/"*** POINTS INPUT • "/I5///5X/ 

2 " MAXIMUM ALLOWED • »/I5//5Xz"*+* RUN ABORTED.") 

1100 F0RMAT13F10.0) 

1110 FORMAT(6I5»8A10) 

1120 FORMAT (1X/BE15.8) 

C ** 

O* 

’ END 

SUBROUTINE CALCORt YERRMX/ ZERRMX/IYLOCJ/IYLOCK/IZLOCJ/IZLOCK) 


79/06/26. 

13*30*37 

C0R01 

309 

C0RD1 

310 

C0R01 

311 

C0RD1 

312 

C0RD1 

313 

C0R01 

314 

COROl 

315 

C0RD1 

316 

COPOl 

317 

COROl 

318 

COROl 

319 

COROl 

320 

COROl 

321 

COROl 

322 

COROl 

323 

COROl 

324 

COROl 

325 

COROl 

326 " 

COROl 

327 

COROl 

328 

OORD1 

329 

C0RD1 

330 

COROl 

331 

C0RD1 

332 

COROl 

333 

COROl 

334 

COROl 

335 

COROl 

336 

COROl 

337 

COROl 

338 

COROl 

339 

COROl 

340 

COROl 

341 

COROl 

342 

COROl 

343 

COROl 

344 * 

COROl 

345 

COROl 

346 

COROl 

347 

COROl 

348 

COPOl 

349 

COROl 

350 

C0RD1 

351 

C0RD1 

352 


*' > 


* ■» 


4 


CREATION RUN 


CAROS ENCOUNTERED IN INPUT UPDATE 1,2-452, 

79/06/26t 

13 *30* 37* # ^ PAGE 

C + + 



C0RD1 

353 

C + + 


THIS SUBROUTINE ITERATES FOR COORDINATE SYSTEM SOLUTION, 

. C0RD1 

354 

C + + 


THE SUBROUTINE IS CALLED FROH COORDC PROGRAM 

CORDl 

355 

C + + 



C0RD1 

356 


COMMON /BXXX/ Y UMXLl, KMX LU i Z { JMXll/ KMXLU , 0( JMXLL) >R (KHXIU/ 

C0RD1 

357 


1 

X{IMXLLJ,YHIMXLL),Z1(IMXLL) 

C0RD1 

358 


COMMON /A XXX/ I LIMIT/ JLIMIT,KLIHIT> IMAX, JMAX> KMAX> IPRT1> IPRT2> 

C0RD1 

359 


1 

IPLTl,IPLT2,ITERMX,Rl,R2,R3,Il,HEDim,HED2<4> 

C0RD1 

360 

C + + 



C0RD1 

361 

C + + 


COMPUTE IN INNER FIELDS 

C0RD1 

362 

C + + 



CORDl 

363 


J1 

• JMAX - I 

CORDl 

364 


K1 

• KMAX - 1 

CORDl 

365 


DO 

200 J • Z,J1 

CORDl 

366 


DO 

200 K • 2, K1 

CORDl 

367 



YPOPO • Y(J,K) 

C0R01 

368 



YP0P1 • YU, K+l) 

C0RD1 

369 



YP0K1 - YCJ,K-1) *"*' 

C0RD1 

370 



YP1P0 • Y (J+1,K) 

C0PD1 

371 



YM1P0 » YU-1, K) 

CORD! 

372 



YM1M1 - YU-1, K-l) 

C0R01 

373 



YM1P1 • YU-1, K+l) 

CQR01 

374 



YP1M1 » Y(J+1,K-1J 

C0PD1 

375 . 



YP1P1 » Y ( J +1, K+l ) 

C0RD1 

376 



ZPOPO « ZU,K) 

C0RD1 

377 



ZP0P1 • Z (J , K+l ) 

C0RD1 

378 



ZP0M1 - Z ( J, K-l ) 

C0RD1 

379 



ZP1P0 » Z(J+1,K) 

C0RD1 

300 ‘ 



,ZM1P0 • Z<J-1,K) 

CORDl 

381 



ZM1M1 ■ Z (J-1,K-1) *. 

C0RD1 

382 



ZM1P1 ■ Z(J-1,K+1> 

CORDl 

383 



ZP1M1 • Z ( J + 1, K-l ) 

CORDl 

384 



ZP1P1 « Z (J+1,K+1) 

CORDl 

385 



QOO • QU) 

CORDl 

386 



RRR • R(K ) 

CORDl 

387 



ALPHA - 0 , 25 ♦ C C YP0P1 - YP0Ml)**2,O ♦ (ZP0P1 - ZP0Hl)**2,0) 

CORDl 

388 •’ 



BETA « 0.25+UYP1P0 - YM1P0) + (YP0P1 - YP0M1) ♦ 

CORDl 

389 


1 

(ZP1P0 - ZM1P0) ♦ (ZP 0 P 1 - ZPOMin 

CORDl 

390 



GAMMA ■ 0.25+ ( ( YP1P0 - YK1P0)**2,0 + (ZP1P0 - ZMIPO)**2.0) 

CORDl 

391 



AJCB1 ■» 0,25* ( ( YP1P0 - YM1P0I * (ZP0P1 - ZP0M1) - 

CORDl 

392 


1 

(YP0P1 - YP0M1) ♦ (ZP1P0 - ZM1P0M 

CORDl 

393 

C++ 



CORDl 

394 . 



YYY1 - <0.5 * ALPHA * (YP1P0 ♦ YM1P0) - 0.25 * BETA * 

CORDl 

395 


X 

( YP1P1 - YP1M1 - YM1P1 ♦ YM1H1J ♦ 0.5 + GAMMA * 

CORDl 

396 



CREATION RUN 


CARDS ENCOUNTERED IN INPUT 


UPDATE 1,2-452, 


CYP0P1 + YPOM1) + 0.25 * AJCBX * AJCB1 *((YP1P0 - YM1P0) * 
QOQ ♦ (YPOPl - YP0M1) *RRR ) ) / (ALPHA ♦ GAMMA) 

ZZZ1 • (0.5 * ALPHA ♦ (ZP1P0 ♦ ZM1P0) - 0.25 * BETA * 

(ZP1P1 - ZP1H1 - ZM1P1 ♦ Z Ml Ml ) ♦ 0.5 ♦ GAMMA * (ZP0P1 

♦ ZP0M1) 4 0.25* AJCB1+AJCB1* (( ZP1PO-ZH1PO) *000 4 
(ZP0P1 - ZP0M1 )*RRR) ) / (ALPHA 4 GAMMA) 

INTERMEDIATE VALUE CALCULATED. COMPUTE NEW VALUE USING 

GAUSS SEIDELL ITERATION 

YYY1 ■ R1 * YYYI 4 (1,0 - Rl) * YPOPO 

ZZZ1 • Rl * ZZZ1 4 (l.o - Rl) * ZPOPO 

FIND MAX ERROR AND ITS LOCATION 

ERR1 • ABS ( YYYI - YPOPO) 

IF (ERR1.LE. YERRMX) GO TO 30 

IYLOCJ • J 

IYLOCK ■ K 
YERRMX « ERR1 

CONTINUE 

ERR1 » ABS ( ZZZ1 - ZPOPO) . - 

IF(ERRl.LE.ZERRMX) GO TO 60 
IZLOCJ ■ J 
IZLOCK ■ K 
ZERRMX ERR1 


60 

.CONTINUE 


C** 



C + + 

SET NEW 

VALUES 

C + + 


Y(J,K) 

* YYYI 


Z( J/K) 

« ZZZI 

C + + 

200 

CONTINUE 



C** APPLY NEUMANN CONDITIONS ON R.H.S. AND UPPER BOUNDARY 

C** 

DO 220 K • 2/K1 

ZUMAX/K) • Z(JMAX-1>K) 

220 CONTINUE 
'C* + 

DO 240 J • 2»J1 
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C0RD1 397 

C0RD1 398 

C0RD1 399 

C0RD1 400 

C0RD1 401 

C0RD1 402 

CORD! 403 

C0RD1 404 

C0RD1 405 

C0RD1 406 

C0RD1 407 

C0RD1 408 

C0RD1 409 

C0RD1 410 

C0RD1 411 

C0RD1 412 " " 

C0RD1 413 

C0RD1 ’ 414 

C0RD1 415 

C0RD1 416 ' “ 

C0R01 417 

C0RD1 418 " ‘ 

C0RD1 419 

C0R01 420 

C0RD1 421 “■ 

C0RD1 422 ‘ " 

C0RD1 423 

C0RD1 424 

C0R01 425 

C0RD1 426 

C0RD1 427 

C0RD1 42 B 

C0RD1 429 

C0RD1 430 

C0RD1 431 

C0RD1 432 

C0RD1 433 

C0RD1 434 ‘ 

C0RD1 435 

CORD! 436 

C0RD1 437 

C0R01 438 ' ' ' 

C0RD1 439 

CORD! 440 


* 


CREATION RUN 


CARDS ENCOUNTERED IN INPUT 


UPDATE 1.2-A52 
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Y ( J.KHAX) ■ YU/KHAX-l) 

2A0 CONTINUE 
C** 

RETURN 

END 

SUBROUTINE IPRTC<ISOLN,KKK>YERRHX>ZERRHX» IYLOCJ>IYLOCK» 

1 IZLOC J> IZLOCK) 

C** 

c** this subroutine prints the coordinate system in y and z 

c** PLANE. ISOLN . 1 INDICATES INITIAL GUESS IS TO BE PRINTED. 

C** ISOLN ■ 2 INDICATES INDICATES CONVERGED OR PARTIAL:. 1 / 

C** CONVERGED SOLUTION IS TO BE PRINTEO. KKK • 0 INDICATES 

C** PARTIALLY CONVERGED SOLUTION AND KKK ■ 1 INDICATES CONVERGED 

C+* SOLUTION 

C + * 

COHHON /BXXX/ Y ( JMXLL.KHXLL ) , Z ( JMXLL/KMXLO* 0( JMXLl ), R (KMXLLJ# 

1 XCIMXLU,Y1(IMXLU>Z1(IMXLL) 

• COHHON / AXXX/ ILIMIT, JL IMIT/ KL I HIT* IMAX/ JHAX* KHAX* I PRT1* IPRT2* 

1 IPLT1*IPLT2*ITERMX»R1*R2*R3*I1»HED1(4)*HE02(A) 

C*+ 

C++ PRINT HEADINGS 

C + * 

IF(ISOLN.EO.l) WRITE ( 6/ 1000 ) 

IF( ISOLN.EQ.2.AND.KKK.EO.O) WRITE(6»1010) 

IF ( ISOLN. EQ.2.AND.KKK.E0.1) WRITE (6# 1020) 

WRITE (6/ 1025 ) HED1 ' ‘ 

WRITE(6>1025) HED2 
WRITE ( 6/ 1030 ) JHAX/KHAX 
IF ( ISOLN. EO.l ) GO TO 50 
WRITE ( 6* 10A0) II 

WRITE ( 6/ 1050) YERRMX* IYLOCJ* IYLOCK 
WRITE ( 6* 1060) ZERRMX,IZLOCJ>IZLOCK 
50 CONTINUE 
C** 

C** PRINT ARRAYS 

c*+ 

DO 100 J • 1>JMAX 
WR 1T6 ( 6f 1070 ) J 

WRITE ( 6/ 1000 ) ((YU/K)iZ(^Kll|K - 1,KMAX) 

100 CONTINUE 

c** 

C++ PRINT J»1 / K*l/ I-l/IHAX 

'C + + 

WRITE ( 6>1200 ) 


C0RD1 

441 

C0RD1 

442 

C0RD1 

443 

C0R01 

444 

C0RD1 

445 

C0R01 

446 

C0RD1 

447 

C0RD1 

448 

C0R01 

449 

C0R01 

450 

C0RD1 

451 

C0R01 

452 

C0R01 

453 

C0PD1 

454 

C0RD1 

455 

C0RD1 

456 

C0RD1 

457 

C0RD1 

458 

C0RD1 

459 

C0R01 

460 

C0R01 

461 

C0RD1 

462 

C0RD1 

463 

C0RD1 

464 

COROl 

465 

C0PD1 

466 

com 

467 

CORDl 

468 

COROl 

469 

CORDl 

470 

CORDl 

471 

CORDl 

472 

CORDl 

473 

CORDl 

474 

CORDl 

475 

CORDl 

476 * 

COROl 

477 

CORDl 

470 

CORDl 

479 

CORDl 

480 

CORDl 

481 

COROl 

482 

CORDl 

483 

CORDl 

484 


CREATION RUN 


C + + 

C+ + 
C++ 
1000 
1010 
1020 
1025 
1030 
10A0 
1050 


CAROS ENCOUNTERED IN INPUT 

WRITE ( 6 , 1210 ) UXm,Yimmiin,I.l,IMAXl 

FORMAT STATEMENTS 


UPDATE 1.Z-A52. 


79/06/26. 13.30.37. 


PACE 12 


FORMAT t 
FORMAT ( 
FORMAT ( 
FDRMATt 
FORMAT ( 
FORMAT ( 
FORMAT ( 


1060 FORMAT { 
1 

1070 FORMATt 
1080 FORMATt 
1 2X 

FORMATt 
FORMATt 

1 


1200 

1210 


C + + 
C + + 
C + + 


C++ 
C + + 
C+ + 
C + + 
C+ + 
C + + 
C + + 
C + + 


C + + 
C + + 
C + + 


»1»//51X, 
"1 ''/MAX* 
"I"/M9X, 
" ",92X» A 
"0",50X,» 
" ",97X,» 
" ",37X, " 
",",I5) 

" "»37X>" 

/IX, »J • 
IX, Ell. 5/ 
, Ell. 5, ZX 
/IX," J-l 
IX, Ell, 5, 
Ell. 5, 


"*** INITIAL SOLUTION ***") 

"*** PARTIALLY CONVERGED SOLUTION *♦*") 

'•**+ CONVERGED SOLUTION » + ♦«) 

A10) 

JMAX • ", I5,AX,*'KMAX ■ ",I5> 

MAX ITERATIONS PERFORMED « «,I5) 

MAX Y-ERROR » ",E15.6,2X,"AT LOC. J,K • ",I5, 

MAX Z-ERROR • ",615.6,2X,«AT LOC. J,K • «',I5> 


" , 1 5 , 


t Y-Z ARRAY)") 


2X,E11.5,2X,E11.5,2X, Ell. 5, 2X, Ell. 5,2 X, Ell ,5, 
,E11.5,2X,E11.5,2X,E11«5) 

, K«l, X*l, IMAX ( X, Y, Z ARRAY)") 

2X,E11.5,2X,E11.5,2X,E11.5,2X,E11.5,2X> 

2X,E11.5,2X,E11.5,2X,E11.5) 


ALL FUNCTIONS COMPLETE 


RETURN 

END 

SUBROUTINE IPLTCtlSOLN, KKK, YERRMX,ZERRHX,I YLOC J,IYLOCK, 

1 IZLOCJ, IZLOCK) 

. TH1S SUBROUTINE IS USEO TO PLOT INITIAL, PARTIALLY CONVERGEO 
OR CONVERGED SOLN. ISOLN • 1 INDICATES INITIAL GUESS 
ISOLN • 2 AND KKK ■ 0 INDICATES PARTIALLY CONVERGED SOLN. 
ISOLN « 2 AND KKK » 1 INDICATES CONVERGEO SOLN. 

GLOBAL COMMON STATEMENTS 

COMMON /B XXX/ Y I JMXLL, KMX LL ) , Z t JMXL L, KMXLL > , 0 ( J HXLL ) , R tKMXLL ) , 

1 X 1 1 HXLL ) , YU IMXLL),Z1(IMXLL) 

COMMON / AXXX / ILIMIT, JLIMIT,KLIMIT, IMAX, JMAX,KMAX,IPRT1, IPRT2, 

1 IPLTl,IPLT2,ITERMX,Rl,R2,R3,Il,HEOim,HED2U) 

LOCAL DIMENSION STATEMENTS 

DIMENSION YDUM t KMXLL ), YOU Ml tJMXLL),ZOUM tKMXLL), ZDUMlt JMXLL ), 
l HE D ( ) 


COROl 

COROl 

CQRD1 

C0RQ1 

CORDl 

CORDl 

C0RD1 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

C0RD1 

COROl 

COROl 

COROl 

COROl 

COROl 

C0RD1 

COROl 

C0RD1 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

COROl 

C0RD1 

COROl 

C0RD1 

CORDl 

COROl 

COROl 

CORDl 

COROl 

COROl 

CORDl 


485 

486 

487 

488 
409 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 
5Q7 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 





CREATION RUN 
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C++ 


CORDl 

529 

C++ 

WRITE HEADINGS. TWO CARDS ARE READ IN THE 

com 

530 

C++ 

MAIN PROGRAM WITH AO COLUMNS INFORMATION ON EACH CARO 

com 

531 

C++ 

WHENEVER IPLT1 OR IPLT2 - 1. 

com 

532 

C++ 


com 

533 


CALL SCRIBE(0.3.1.5.0.15.0,5.HED1.90,.A0.9) 

CORDl 

534 


CALL SCRIBE(0.55.1.5.0«15.0.5.HED2.90*.A0.9) 

com 

535 

C + + 


CORDl 

536 

C+ + 

PRINT HEADINGS FIRST 

CORDl 

537 

c** 


CORDl 

538 


IF ( ISOLN. E0.2 > 60 TO 50 

CORDl 

539 


HED(l) • IOH INITI A 

C0R01 

540 


HE0(2) • 10HL GUESS — 

CORDl 

541 


HED C 3 ) « 1H- 

CORDl 

542 


CALL SCRIBE(0. 70/3,0/ 0,1/ 0.5. HED. 90 *0.21. 9} 

CORDl 

543 


GO TO 100 

CORDl 

544 

C + + 


CORDl 

545 

C + + 

PARTIALLY CONVERGED OR CONVERGED SOLN 

CORDl 

546 

C+ + 


CORDl 

547 

50 

CONTINUE 

CORDl 

548 


IF(KKK.EQ.l) GO TO 75 

CORDl 

549 


HED ( 1 ) • IOH PARTIA 

CORDl 

550 


HE D ( 2 ) » 10HLLY CONVER 

CORDl 

551 


HED ( 3 )- ■ 10HGED SOLUTI 

CORDl 

552 


HED(A) » 7H0N 

CORDl 

553 


CALL SCRIBE (0,70. 3. 0.0. 1.0. 5. HED. 90.0.37.9) 

CORDl 

554 


GO TO 100 

CORDl 

555 

75 

CONTINUE 

CORDl 

556 


.HED(l) » IOH CONVER 

CORDl 

557 


. HED ( 2 ) - 10HGED SOLUTI 

CORDl 

558 


H ED C 3 ) ■ 6H0N 

CORDl 

559 


CALL SCRIBE (0,70/ 3. 0.0, 1.0. 5. HED. 90. 0.26. 9) 

CORDl 

560 

100 

CONTINUE 

CORDl 

561 

C + + 


CORDl 

562 

C + + 

WRITE ARRAY SIZES 

CORDl 

563 

C + + 


CORDl 

564 


ENCODE ( AO. 1000. HED ) JMAX.KMAX 

CORDl 

565 


CALL SCRIBEtO, 05/3,0.0.1.0, 5. HE 0.90. 0.27. 9) 

CORDl 

566 

C+ + 

. 

CORDl 

567 

C + + 

CALL NFRAME TO ADVANCE FRAME AFTER WRITING HEAOINGS 

CORDl 

568 

c** 


CORDl 

569 


CALL NFRAME 

CORDl 

570 

'C + + 


CORDl 

571 

C + + 

FIND MIN MAX VALUE FOR PLOTTING 

CORDl 

572 


RE ATI ON RUN 


CAROS ENCOUNTERED IN INPUT 


. UPDATE 1.2-452 


C + + . 

YMIN • 10000.0 

ZMIN • 10000.0 
YMAX ■ -10000.0 
ZMAX - -10000.0 
DO 150 J » l/JMAX 
DO 150 K « l/KMAX 

IF(Y(J/K).LT .YMIN) YMIN • Yt J/K) 

IFCY( J,K) .GT.YMAX) YMAX ■ Y(J»K) 

IF(ZtU,K) .IT. ZMIN) ZMIN » Z(J,K> 

IF(Z(J,K) .GT.ZMAX) ZMAX • Z(J,K) 

150 CONTINUE 
C++ 

C + + PLOT DATA. ZETA LINES FIRST 

C + + 

DO 210 J « l/JMAX 
DO 200 K ■ l/KMAX 
YDUM(K) • Y(J,K) 

ZOUM(K) ■ Z( J/K) 

200 CONTINUE 

CALL INFQPLT(C,KMAX,YDUM,1,ZDUM,1, YMIN/ YMAX, ZHIN/ZMAX, 

1 1 .O/O/HEO/O/HEO/O/ 5.0/ 5 .0/0.0/ 2.0) 

210 CONTINUE 
C + + 

C++ PLOT ETA LINES 

C + + 

KK • 0 

DO 250 K * 1/KMAX 

IF(K.EO.KMAX) KK » 1 
. DO 240 J ■ l/JMAX 

YDUM1IJ) » Y(J,K) 

ZDUhl(J) ■ ZIJ/K) 

240 CONTINUE 

CALL INF0PLT(KK,JMAX,YDUM1,1,ZDUM1/1, YMIN, YMAX, ZMIN, ZMAX, 
1 1.0,0,HED,0,HED, 0,5. 0,5. 0,0. 0,2.0) 

250 CONTINUE 
C++ 

RETURN 

C + + 

1000 F0RMATI7HJMAX • , IS, 10H KKAX « ,15) 

C** 

END 


A 
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C0R01 

573 

CORDl 

574 

CORD! 

575 

CORDl 

576 

CORDl 

577 

CORDl 

578 

CCRD1 

579 

CORDl 

580 

CORDl 

581 

CORDl 

582 

CORDl 

583 

CORDl 

584 

CORDl 

585 

CORDl 

586 

CORDl 

587 

CORDl 

588 

CORDl 

589 

CORDl 

590 

CORDl 

591 

CORDl 

592 

CORDl 

593 

CORDl 

594 

CORDl 

595 

CORDl 

596 

CORDl 

597 

CORDl 

598 

CORDl 

599 

CORDl 

600 

CORDl 

601 

CORDl 

602 

CORDl 

603 

CORDl 

604 

CORDl 

605 

CORDl 

606 

CORDl 

607 

CORDl 

608 

CORDl 

609 

CORDl 

610 

CORDl 

611 

CORDl 

612 

CORDl 

613 

COPPi 




B. Job Control Cards 


A list of job control cards needed to run this pro- 
gram are included in the comment section of the main program 
COORDC. 

C. Input Explanation 

The instructions for input preparation are included 

t 

in the comment section of the main program COORDC. 

D. Sample Input 

Sample input listing follows on the next page. 

E. Sample Output 

The sample output generated by the sample input of 
section D follows sample input. 
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PROGRAM NO- 

CODED BY 

DIVISION 


LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


DATE- 

PAGE- 


-SECTION- 


JOB ORDER. 



NASA- Langley Form 67 (MAR 69) 


NOTE. WRITE NUMBERS 10, LETTERS lj/UG7C, SYMBOLS / 


*US GOVLRNMEI 


1 











NASA-Lang.ley Form 67 (MAR 69) 


NOTE: WRITE NUMBERS 10, LETTERS I $ U G -Z C, SYMBOLS / . , 


*US UOVEl^ 






PROGRAM NO- 

CODED BY 

DIVISION 


LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


PAGE- 


-SECTION. 


JOB ORDER. 


STATEMENT 

NUMBER 


1 2 3 4 5 


L 0 ( . 0 

0 . 0 

. L L ... 

,_P, . o, 

0.0 

, 0 , . 0 
, 0 , . 0 

0 . 0 

L. .1 

0 . 0 


.i i-i 9 l.: 

i i 


...iPu 


, . , °, 



NASA-Lungley Form 67 (MAR 69) 


NOTE: WRITE NUMBERS 10, LETTERS Ij/UG7C, SYMBOLS / . , 


*U b GOVE. H N ML 






PROGRAM NO.- 

CODED BY 

DIVISION 


LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


-SECTION- 


JOB ORDER 


STATEMENT 

NUMBER 


1 2 3 4 5 

. Q, • 

.!.... L 

. , , 9l- 

L 

, i C, 0 0 


-CONTINUATION 


FORTRAN STATEMENT 


L t ! ? IP. !!J2 13 J4 J5 16 .1_7_J8 J9 20.21 .22 2?„24 2£ ..24_2I 20 22_10 J4-35 J6-J7. 40 _ 41.42 43 44-4i-46_42_ 

Qi_Zi_6u0j |.Q ua_P -0-i IlI 1 3 l8_i — i-Qj-a.-i.Q- > 3 1 0 * 1 ■ 1 1 3 > S_i — i — ». — » — i — i — i — i — i — » — i — 4 — * — i — j — 1 — 1 — u 

Zj-jLiPLiPi^PxllP.i 1 . 3 , 8 I ., 0 , . , 0 ,310,1, 1 , 3 , 8 


J... L_ _L- . 

. .L. _i„l 

. .1 1 .. J ..... 

I .... L — L. .. 

. _ l l l — 

l 1 _.a _ 

_. L-..1 l. 


7,_6, Oj 

o 

o 

CO 


A o,o| 

0 l’ j °i 3 

0,1, 1,3, 8, 

8^4, 0{_ 

0, .,0,3 

0,1, 1,3, 8, 

e r If 

; i l | 

M, W , 

W.I.T.H, | F 

I.N.Aj Tl 

e l .sjr 

S , T, E, Mi .1 


X- _ 1 

J„ .1 . .. 

J _ . 1 . - 
.. I .J .. 

L 

i L . 


j 1 l I 

\ _ 1 J L 

i — i -l—i_ 

j J .1 1 

j I 

J A„_J J 1 


. 1 t 1 1 . i ii. I 

Si T, E, M, * 1 , S , I , Z , Ei - , , 3, Oi Xi 2 . 3 ,X,2 ,3i 

I 1 __i I L L 1 1 4 1 I 1 J 1 .■■■■.] 1 I I 1 L 

„j 1 1 ,„i ...1 1 1 1 1 L 1 1 1 — 1 — I 1 1 1 1 L 

1 1 j 1 1 t... . 1 l__j L 1 1 j 1 I 1 1 1 1 L_ 

X--A l 1 1 I -L_ J J . 1 , 1 . ..I U I l I I L. ..J 1 _ 

l L_J i L 1 l J I 1 J. — J J A [ I L I l 1 _ 


L- -1... ... l_ 1. i !. I i 1. I 1 L- 

J t J ! A J l J l J I — 

. J i__J ! J J. , L L I 

1 1 1 ± I > 1 > | I i 1 

1 L -J 1 1 1 1 L I L L. 


NASA- Langley Fom G'J (MAR 69) 


NOTE: WRITE NUMBERS 10, LETTERS irfUGiC, SYMBOLS / . , * 






E . SAMPLE OUTPUT 


+** INITIAL SOLUTION +♦♦ 
CORNEP FLOW WITH FILLET 
COORDINATE SYSTEM* SIZE»30X23X23 


JMAX • 23 KMAX ■ 23 


J • 1 

.30119E-01 
• 91 3335-02 

0. 

0* 

(Y-Z ARRAY) 
.301192-01 
•73335E-01 
.218922+00 
.999052+00 
•67209E+00 

.267232-01 
. 93390E-03 

0. 

0. 

0. 

.33673 E-01 
• 93379 E-01 
.259335+00 
•99200E+00 
.712095+00 

.229175-01 

0. 

0. 

0, 

0. 

.387325-01 
.118075+00 
.303035+00 
• 53930E+00 
.750005+00 

• 16805E-01 

0. 

0. 

0. 

.969875-01 
.197392+00 
• 398865+C0 
.585382+00 

.103705-01 

0. 

0. 

0. 

.578205-01 
.180925+00 
.396105+00 
* 62979E+00 

J * 2 

.336735-01 
* 7 0219E-O2 
•370G5E-02 
« 37035E-02 
. 37085 E-02 

(Y-Z ARRAY) 
•267232-01 
•701972-01 
.215922+00 
.992612+00 
. 67168 E + 00 

.30299E-01 
.919032-02 
• 370G5E-02 
.370655-02 
♦ 3706 5E-02 

.302995-01 

.902865-01 

• 25701 E + 00 

• 99079 E + 00 
.711925+00 

.260195-01 
.370855-02 
.370855-02 
.370055-02 
• 37085 E-02 

.353825-01 

.115105+00 

.300925+00 

•53831E+00 

•75000E+00 

.209305-01 
.370855-02 
• 37085E-02 
. 3708 5E-02 

.931735-01 

.199505+00 

•39697E+00 

.589605+00 

•190285-01 
.370855-02 
•37085E-02 
. 370855-02 

.595595-01 

•178295+00 

‘.399935+00 

.629185+00 

J • 3 

.38732E-01 

• 13Q63E-0 1 

• 697Q8E-02 
.897885-02 

• 8.978 bc-02 

(Y-Z ARRAY) 
• 22917 E-01 
. 66100E-01 
•21279E+00 
•99078E+00 
•6712 12+00 

.353622-01 

.990755-02 

.897665-02 

•89788E-02 

.897686-02 

.260195-01 
.863585-01 
.259085+00 
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.750005+00 

. 431 25E + 00 

.750005+00 

*481216+00 

•75000E+00 

• 53049E+00 

♦75000E+00 

« 57849E + 00 

•75000E+00 

.624716+00 

• 75 000E + 00 

.668785+00 

.750005+00 

.71051 E+00 

.750006+00 

•75000E+00 





J«1>K*1,X»1 

/ 1MAX ( X> Y> Z 

ARRAY) 








0. 

.301145-01 

. 301 14 E-01 

• 68190 E-02 

.301146-01 

•30114E-01 

•13640E-01 

• 30114E-01 

• 30114E-01 

* 

♦204805-01 

. 30114E— 01 

. 30 1 14E-01 

.273706-01 

.301146-01 

.301146-01 

.344606-01 

• 30114E-01 

•30114E-01 


• 4 212 0E-01 

.301146-01 

• 30 1 1 4E-01 

* 5 1190E-01 

.301146-01 

.301146-01 

•63230E-01 

• 30114E-01 

• 30114E-01 


* 8Q520E-01 

.301145-01 

. 30 1 14 E-01 

• 10560 E + 00 

.301146-01 

« 30 114E-01 

♦ 14120E+00 

• 30114E-01 

. 30 114E-01 


.176606+00 

•30114E-01 

.301146-01 

.212006+00 

♦30114E-01 

.301146-01 

. 24740E +00 

. 30114E-01 

• 30114E-01 


.282805+00 

.301145-01 

. 30 1 14 E-01 

* 31820 E + 00 

. 301 14 E-01 

* 30 114 E-01 

• 35360E+00 

.30114E-01 

♦ 30 114E-01 


.33900E+00 

.301146-01 

.301145-01 

• 4 24 4 OE + 00 

. 301 1 4E-01 

♦30114E-01 

.459606+00 

• 301 1 4E-01 

•30114E-01 


.495205+00 

.301145-01 

.301145-01 

• 53060E +00 

•30114E-01 

•30114E-01 

♦ 56600E +00 

.30114 E-01 

.301145-01 


.601405+00 

.301145-01 

♦301145-01 

.636805+00 

.301145-01 

♦30114E-01 

.672205+00 

•30114F-01 

. 30114E-01 


• • 7Q760E +00 

. 301 14 E-01 

.301145-01 

♦743006+00 

.301145-01 

.301145-01 

•77040E+00 

• 30114E— 01 

.30114E-01 


ITER NO 

IYLOCJ* IYIOCKj YERRKX* IZIOC J/IZ tOCK» ZERRKX • 50 

13 9 

.7161946-04 

9 13 

•715067E-04 




*** PARTIALLY CONVERGED SOLUTION *+* 

CORNER FLOW WITH FILLET 
COORDINATE SYSTEM# SI7E-30X23X23 

JMAX « 23 KMAX » 23 

MAX ITERATIONS PERFORMED • 50 

MAX Y-ERRQP * .71619AE-0A AT LOC. J/K • 13/ 9 

MAX Z-ERROR • .715067E-0A AT LOC# J>K • 9, 13 


J ■ 1 

• 301 1 A 6-01 

(Y-Z ARRAY) 
• 3011 A6-01 

• 26723E-01 

#336736-01 

•22A17E-01 

.A1333E-02 

.733356-01 

.A3390E-03 

#93379 E-01 

0. 

C. 

•218A2E+00 

0. 

.259336+00 

0. . 

0. 

» AAA05E + 00 

0. 

♦A9200E+00 

0# 

0. 

• 6720A 6 + 00 

0. 

.712096+00 

0. 

J • 2 

.336736-01 

(Y-Z ARRAY) 
* 26723 E-01 

• 31702E-01 

•31701 E-01 

•27769E-01 

.631996-02 

. 7272 5E-01 

» A2756E-02 

.927A2E-01 

.360A3E-02 

.356166-02 

• 217 6 A E +00 

• 35 6 1 2E-02 

• 2 588 1 E + 00 

• 3 5609E-02 

.356076-02 

. AA375E+00 

.356076-02 

. A9175 E + 00 

♦35607E-02 

.356076-02 

.671986+00 

.356076-02 

.71206E+C0 

.356076-02 

J * 3 

.367325-01 

(Y-Z ARRAY) 
• 22A 17 6-01 

♦376176-01 

•27767E-01 

♦3A7fc36-01 

.1A5U7E-01 

* 717 AO E-01 

.905916-02 

.91793 E-01 

.877866-02 

.662756-02 

. 2 1 705 £ +00 

• 86 2 A9E-02 

.258096+00 

.062376-02 

.862276-02 

♦A A 3356 +00 

. 86227E-02 

# A91 A 2 E + 00 

•86227E-02 

.362266-02 

.671896+00 

•86226E-02 

.712026+00 

. 86228E-02 

J * A 

, A6A 8 7 E-Q 1 

(Y-Z ARRAY) 
.168056-01 

* A5 7 A 16-01 

• 21862 E-01 

. A39A 56-01 

« 2 A717E-01 

• 697 36E-01 

.189306-01 

.901 A3E-01 

•169A8E-01 

•16A15E-01 

.215876+00 

• 16 A02E-01 

♦257056+00 

. 16396E-01 

.163926-01 

. AA276E+00 

.163916-01 

•A909AE+00 

• 16391E-01 

.163925-01 

.671776+00 

.163926-01 

• 7 1197 E + 00 

• 16392 E-01 

0 • 5 

•07e20£-01 

(Y-Z ARRAY) 
.103706-01 

.571876-01 

•1A951E-01 

.559876-01 

.398766-01 

• 6570A6-01 

.329936-01 

.871766-01 

. 29A90E-01 

.278706-01 

•21A22E+00 

♦ 27 818 E-01 

.255596+00 

• 2779A E-01 

• 2777 A £-0 1 

.AA1975+00 

, 2777 A6-01 

.A90296+00 

,277736-01 

,277755-01 

,671606+00 

#277756-01 

.711696+00 

.277756-01 

J • 6 

,733366-01 

(Y-Z APRAY) 
, A1 333 6-02 

•727266-01 

•83196E-02 

♦717A2E-01 

.593586-01 

.593356-01 

.523696-01 

.023526-01 

•A7A296-01 


#387326-01 

• 11807E + 00 
.303036+00 

# 53930E + 00 
#750006+00 

•16805E-01 

0# 

0# 

o. 

# A6A87E-01 
•1A73AE+00 
•3A886E+00 
• 50538E+00 

#103706-01 

0# 

0# 

0# 

.578206-01 
#180926+00 
•396106+00 
# 6297A 6+00 

#376166-01 
* 11739E+00 
•30256F+00 
#539116+00 
•7500CE+00 

♦ 2188AE-01 
•35675E-02 
#356006-02 
•35607E-02 

• A57A0E-01 
• 1 A666F+00 
•3A8A56+00 
#585236+00 

• 1A952E-01 
#356316-02 • 
• 35607E-02 
. • 35607E-02 

#571866-01 

• 18028 E+00 

• 39575 E + 00 
•6296AE+00 

#3 A759E-01 
•116A2E+00 

• 30193E + 00 

• 5308 AE + OO 

• 7 5000E + 00 

•29291 E-01 
•8653AE-02 
• 86231E-02 
• 86227E-02 

• A39A0E-01 
♦1 A572E+00 

• 3A790E + 00 
.505036+00 

• 21B59E-01 
• 86337E-02 
• 86228E-02 
• 86228E-02 

• 55983E-01 

• 179A1E + 00 

• 39527 E+00 

• 62950E+00 

.292856-01 
.11 A90E+00 
.301016+00 
♦ 53 8A6E + 00 
.750006+00 

#397516-01 
. 16537E-01 
• 1639AE-01 
#163926-01 

• 397A1E-01 
• 1AA30E + 00 
.3A7096+00 
•58A7A6+00 

• 32753E-01 
• 16AA5E-01 
• 16392E-01 
.163926-01 

• 531356-01 
.170106+00 
•39A586+00 
#629296+00 

•21D5A6-01 

.112506+00 

.2997AE+00 

,537956+00 

,750006+00 

• 531 A56-01 
• 283A06-01 
.277836-01 
♦2777AE-01 

# 327A2 6-01 
#1 A22CE+00 

• 3A 599E + 00 
•58A35E+00 

• A7573E-01 
# 27993E-01 
, 27777E-01 
•2777AE-01 

# A7556E-01 
#176256+00 

• 3936A E + 00 
•62902E+00 

•1 A506E-01 
• 10 885E+00 

. 697AAE-01 
• A 50566-01 

#2 A7106-01 
.139266+00 

#657216-01 

.cun 

,398626-01 

1 7l74C*ftn 



I 


4 






. 




... 


.437336-01 

• 21205 5 + 00 

* 43 567E-01 

• 25370E+00 

• 43490E-01 

♦298106+00 

•43452E-01 

.344585+00 

.434335-01 

• 43425 E-0 1 

.440975+00 

• 43421 E-01 

•48947E+00 

.43420E-01 

•53730E+00 

♦43421E-01 

.583876+00 

.434235-01 

,4342*5-01 

.671395+00 

.434255-01 

•71179E+00 

•43425E-01 

♦75000E+00 




J * 7 

{Y-Z ARRAY) 








.933795-01 

.433905-03 

•9274?E-01 

• 4 27 56E-02 

•91794E-01 

.985075-02 

.901486-01 

.1 8928E-01 

. 0719OE-O1 

•823735-01 

.523505-01 

•76281E-01 

• 762 55E-01 

•70797E-01 

.103995+00 

•67247E-01 

.135406+00 

• 65 414E-01 

•645475-01 

.209375+00 

•64137E-01 

.251415+00 

.639365-01 

.296125+00 

.638365-01 

• 342P9E+00 

.637875-01'* 

.637635-01 

. 43979 E + 00 

.637535-01 

• 48852E + 00 

.637505-01 

.536565+00 

.637515-01 

• 58332E + 00 

. 63755E-01 . 

.637585-01 

.671165+00 

.637605-01 

•711695+00 * 

*637605-01 

,750005+00 




J * e 

(Y-Z ARRAY) 








.116075+00 

0. 

• 11739E+00 

.360435-02 

•11642E+00 

.877855-02 

♦114905+00 

.169476-01 

•11251E+00 

• 1C 3675 + 00 

. 47416E-01 

.104015+00 

.707755-01 

,988935-01 ' 

•98866E-01 

♦947635-01 

.131035+00 

.920965-01 

.905955-01 

.206305+00 

.697995-01 

,248795+00 

•89384E-01 

• 29 390E + 00 * 

f89 1686-01 

•34100E+00 

, 89058 E-01 

• 8900 <• £-01 

,438485+00 

• 88979 E-01 

, 48746 E + OO 

•68972E-01 

•535745+00 

.889735-01 

•582716+00 

.889785-01 

♦889846-01 

.670905+00 

.889895-01 

.711575+00 

•88989E-01 

.750005+00 




J • 9 

(Y-Z ARRAY) 








. 1 47345+00 

0. 

♦14666E+00 

• 3 5674 E-02 

* 14 572E + 00 

.865325-02 

.144305+00 

.165366-01 

.142216+00 

.139275+00 

♦450495-01 

• 13 542E+00 

,672325-01 

,131055+00 

.947415-01 

.126985+00 

.126966+00 

.123865+00 • 

.121795+00 

,203095+00 

. 12056E+00 

.246045+00 

• 11986E + 00 

•29154E+00 

• 11947E+00 

* 33900E + 00 

. 11926E + 00 

.119165+00 

• 43711 E + QO 

.119115+00 

,486375+00 

.119095+00 

• 53489 E + OO 

• 11909E+00 

• 58209E+00 

•11910E+00 

• 11911E+00 

.670645+00 

•119115+00 

• 711 4 5E + 00 

•11911E+00 

•75000E+00 




J » 10 

(Y-Z ARRAY) 








.180925+00 

0. 

. 1802 8 E + OO 

. .356315-02 

♦179415+00 

.863365-02 

.178105+00 

.164446-01 

♦176246+00 

.1737 45 + 00 

.441075-01 

•170585+00 

•65404E-01 

,166945+00 

.920795-01 

.163275+00 

• 1238 45+00 

* 16009E+00 

.157715+00 

. 20012 E + OO 

• 15 61 15 +C0 

• 243 39E + 00 

# 15512E + 00 

• 20924E + 00 * 

.154545+00 

♦337046+00 

• 15421E+00 

.154045+00 

• 43576E +00 

.153955+00 

.485305+00 

.153925+00 

.534075+00 

•15391 E+OO 

.581486+00 

♦153925+00 

•15393E+00 

. 67039 E +00 

.153945+00 

,711345+00 

,153945+00 

♦75000E+00 




J * 11 

(Y-Z ARRAY) 








*218425+00 

0. 

• 21784 E +00 

.356175-02 

.217045+00 

•86274E-02 

•21587E+00 

.16414E-0X 

.214226+00 

♦ 21 205E + QQ 

• 437 30E-01 

.209385+00 

.645415-01 

.206315+00 

.905835-01 

.203115+00 

.1217eE+00 

.200146+00 

.197715+00 

.197695+00 

.195925+00 

,241005+00 

.194715+00 

.287175+00 

•193965+00 

.335245+00 

.193515+00 

•19326E+00 

• 4 34 52 E + 00 

• 19 313E +00 

.484315+00 

♦ 19307E+00 

.533315+00 

.193065+00 

. 58093E+00 

* 19307E+00 * 

• 1930 EE + 00 

•670165+00 

*.193095+00 

.711235+00 

•19309E+00 

.750005+00 




J ’ 12 

(Y-Z ARRAY) 








,259336+00 

0. 

• 25 680 £ +Q0 

•356125-02 

• 25809E +00 

.862495-02 

•25704E+00 

• 16402E-01 

.25559E+00 

* Z537OE + 0O 

# ^ 3566 E— 01 

• 25H0E+C0 

#641 33 £-01 

•2480OE+OO 

1097916-01 

•246046+00 

•120556+00 

•243406+00 

f 24UQE + 00 

•195905+00 

•239285+00 

*239 27E + 00 

•237985+00 

.285466+00 

• 23711E+00 

♦333726+00 /’ 

.236575+00 

.236265+00 

,433445+00 

• 23 609 E +00 

•483456+00 

•236016+00 

.532656+00 

.235995+00 

♦580456+00 p 

'.235996+00 

.236006+00 

.669976+00 

,236016+00 

,711146+00 

♦236016+00 

#750006+00 


. .. 


J ■ 13 

(Y-Z ARRAY) 






• 


,303035+00 

0* 

•30256E+00 

*356095-02 

.301935+00 

.862375-02 

.301016+00 

.163966-01 

.299736+00 

,293095+00 

.43489E-01 

♦296125+00 

*639345-01 

. *293896+00 

,893795-01 

. .291546+00 

tll98«+00 — 

# ?flO 94C*AA 


« X f 3 T*lt»VU 


.392436+00 

•628696+00 


•329836-01 
•170576+00 
• 39101E+00 
•62830 E +00 


.294865-01 

•16692E+00 

*389426+00 

•62788E+00 


.283375-01 
•163245+00 
♦ 38776 6 + QO 
•62746E+00 


•27991E-01 
•160075+00 
• 386UE + 00 
.627045+00 


.278696-01 
•1 5769F+00 
•384596+00 
•62667E+00 


• 270.176-01 

• 156106+00 

• 38329 E+OO 
.626346+00 


.277946-01 

. m \ RUl . CaAA 





• 28717E+00 

• 199706+00 

.285576+00 

•237976+00 

.285186+00 

.282356+00 

. 53257E+00 

, 28215E+00 

•48275E+00 

.282056+00 

.282026+00 

•66981E+00 

.282026+00 

•711076+00 

.282026+00 

J » 19 

.358866+00 

(Y-Z ARRAY) 
0, 

.35855E+00 

.356086-02 

.357896+00 

.355576+00 

•53551E-01 

.35288E+00 

•638356-01 

•35099E+00 

» 335256+00 

. 19395E+00 

.333736+00 

.237116+00 

.332556+00 

• 330736+00 

, 53191 E+00 

.33052E+00 

• 48222 E+00 

•33051E+00 

• 33036E+00 

.66969E+00 

» 33036E+00 

•71101 E+00 

•330366+00 

J • 15 

» 396106+00 

(Y-Z ARRAY) 

0. 

.395756+00 

•356C7E-02 

.395276+00 

♦39253E+00 

,535332-01 

.391016+00 

.637866-01 

.389526+00 

.36960E+00 

, 19350E+00 

.38329E+00 

♦ 23657E + 00 

.382256+00 

.360606+00 

« 531 55 E+00 

. 33039 E +00 

.481346+00 

.380286+00 

•380216+00 

« 66960E + 00 

, 38021E+00 

.71097E+00 

.360216+00 

J « 15 

.599056+00 

(Y-Z ARRAY) 

0. 

.553756+00 

• 35607 E-02 

.55335E+00 

.99096E+00 

. 53525E-01 

•■53978E+00 

♦ 637 63E-01 

.53857E+00 

• 939526+00 

. 19326 E + OO 

• 533 55E + 00 

♦23626E+00 

.532576+00 

.‘•3115E + 00 

.531156+00 

.53096E+00 

•48158E+00 

•530P5E+00 

• 93077E+00 

.669596+00 

. 53077c +00 

• 71095 E + 00 

.53077E+00 

J * 17 

.99200E+00 

(Y-Z ARRAY) 
0. 

.591756+00 

• 35607 E-02 

. 59152 E + 00 

.989976+00 

.5 35 2 IE— 01 

. 5085 1E + 00 

♦ 63753E-01 

.58756E+00 

.989316+00 

, 19313E+00 

.58355E+00 

♦23609E+00 

•50275E+00 

.561596+00 

, 53096E+00 

.50153E+CO 

•40142E+C0 

.581336+00 

»9bl 25E+00 

.669506+00 

.58125E+00 

.710936+00 

.58125E+00 

J * 18 

. 53930E + 00 

(Y-Z ARRAY) 

0. 

.53911E+00 

• 3 5607 Er02 

.53885E+00 

.53730E+00 

.53521E-01 

•53665E+00 

• 63750 E-01 

•53573E+00 

•53331E+00 

. 1 9307c +00 

.53265E+00 

• 23601 E + 00 

•53212E+00 

•63121E+CO 

,530856+00 

•5310BE+00 

♦ 4 81 33E + 00 

.531016+00 

• 53093E +00 

. 6695 7E + 00 

, 53093E +00 

♦71092E+00 

. 53093E + 00 

J ■ 19 

(Y-Z ARRAY) 




. 585306+00 

0. 

.585236+00 

.356076-02 

.585036+00 

,583876+00 

.53521E-01 

.583316+00 

•637526-01 

.582716+00 

.58093E+00 

, 19306E+00 

. 58055E+00 

• 23599 E + 00 

.580066+00 

• 579906 +00 

. 53080E+00 

. 579306+00 

• 4 01 2 8 E + 00 

.57925 E+00 

.579186+00 

.66956E+00 

,579176+00 

♦710916+00 

.579176+00 

J « 20 

,629756+00 

(Y-Z ARRAY) 

0. 

.629656+00 

.356076-02 

.62950E+00 

,626696+00 

.939236-01 

.620306+00 

•637596-01 

.627886+00 


.284186+00 

.532126+00 

•75000E+00 

*283286+00 

.232016+00 

*332536+00 

•580066+00 

.282706+00 

.282016+00 

.382256+00 

«62b086+00 

•86231E-02 

•89166E-01 

♦28328E+00 

•53171E+00 

•75000E+00 

.357096+00 

.339006+00 

•331686+00 

•33036E+00 

.163956-01 

.119576+00 

.331676+00 

.579766+00 

.355986+00 
. .337056+00 
.331096+00 
.330356+00 

.277826-01 

.155536+00 

•381586+00 

.625806+00 

•862286-02 
•890576-01 
• 28269E+00 
♦53141E+00 
•75000E+00 

.395586+00 

.387756+00 

.381586+00 

.380236+00 

.163926-01 

.119266+00 

.331096+00 

.579556+00 

-.393636+00 

.386116+00 

.380956+00 

.380216+00 

.277776-01 

.155216+00 

.380956+00 

•62573E+00 

•86227E-02 
• 89003E-01 

• 28 234E + 00 

• 5312 1E + 00 
•75000E+00 

.552766+00 
. 53710E+00 
.531916+00 
.530306+00 

. 16392E-01 
.119166+00 
.330736+00 
.579506+00 

.551976+00 
.535766+00 • 
.531556+00 
.530776+00 

.277756-01 

.155056+00 

.380606+00 

.62563E+00 

• 8 6227E-02 
♦88980E-01 

• 28 214E + 00 
•53108E+00 

• 7500CE + 00 

.590956+00 

.586376+00 

.582226+00 

.501286+00 

.163916-01 

.119116+00 

.330526+00 

.579306+00 

.590296+00 

.585296+00 

.581856+00 

.581266+00 

.277756-01 
*15395 F+00 
•380396+00 
.625566+00 

•86227E-02 
.869726-01 
• 28205E+00 
« 53101E + 00 
•75C00E+00 

,530566+00 
•53569E+00 
. 53171E+00 
. 53096E+00 

•16391F-01 

.119096+00 

.330516+00 

.579256+00 

.537956+00 
.535066+00 
.531516+00 
, 53095E+00 

.27775E-01 

.153926+00 

.380286+00 

.625526+00 

.962276-02 
.889736-01 
• 2 8 201E + 00 
.530966+00 
♦750006+00 

.585756+00 

.582086+00 

,579766+00 

.579216+00 

.163926-01 
.119096+00 
.330366+00 
.579216+00 . 

.5P535E+09 
.581506+00 
, 57955E+00 
.579196+00 

.277756-01 

.153916+00 

.380236+00 

.62550F+00 

.862286-02 

.889766-01 

.62929E+00 

.627556+00 

.163926-01 

.119106+00 

.629026+00 
•677ftAC+nn • 

.277756-01 
-1 Mc?e*nr> 



- - - . .. -w- - tvttvtCtUU 


.626665+00 

• 19307E+00 

•6263AE+00 

♦23599E+00 

*626085+00 

.282015+00 

•62588E+00 

•33035E+00 

.625735+00 

•625A8E+00 

.62563E+Q0 

. A3077 E+00 

*625566+00 

• A8126E+00 

.625525+00 

•5309AE+00 

•62550E+00 

•57919E+00 

• 625A7E+00 

* 669A5 E + 00 

♦ 62 5 A7E+00 

•71090E+00 

•625A7E+00 

•75000E+00 



J • 21 

(Y-Z ARRAY) 








• 6720 AE + OO 

0* 

* 67198E+00 

•35607E-02 

•671896+00 

•86228E-02 

•67177E+00 

•16392E-01 

.671605+00 

• 67139E +00 

.A3A2AE-01 

.671165+00 

• 63758 E-01 

•670905+00 

.869655-01 

•6706AE+00 

.119116+00 

•67039E+00 

.67016E+00 

• 1 9306E + 00 

• 66997E + 00 

.236006+00 

•66981E+00 

.282025+00 

•66969E+00 

.330366+00 

.669605+00 

« 6695 AE+OO 

* A3077E +00 

* 669 50 E >00 

. A012 5 E + 00 

* 669A7E + 00 

.530936+00 

•669A6E+00 

.579165+00 

. 669A5E+00 

« 669 A Ac+OO 

• 669AAE+00 

. 669 A AE+OO 

•71090E+00 

•669AAE+00 

•75000E+00 



J ■ 22 

(Y-Z ARRAY) 








* 71 209 E +00 

0. 

*712066+00 

.356076-02 

•71202E+00 

•86226E-02 

•71197E+00 

•16392E-01 

. 71189E+00 

•71179E+00 

.A3A25E-01 

.711696+00 

♦63760E-01 

.711575+00 

.38989 E-01 

•711A5E+00 

.119115+00 

.7113AE+00 

.711235+00 

• 19309 E + 00 

. 711 1 A£+00 

• 23601 E + 00 

•71107E+00 

•20203E+00 

•71101E+00 

.330365+00 

.710975+00 

♦71095E+0Q 

♦ A3077E+00 

.710936+00 

•A81 256+00 

•710925+00 

• 53093E + 00 

♦71091 E+OC 

.579176+00 

•71090E+00 

♦71090E+00 

.669 A A E +00 

•71090E+00 

.710906+00 

.710905+00 

•75000E+00 



J ■ 23 

(Y-Z ARRAY) 








. 75000E+00 

o. 

• 75000E +00 

•35607E-02 

.750005+00 

• 8 6228E-02 

•75000E+00 

•16392E-01 

.750005+00 

.75 OOCE + OO 

* A3A.25E-01 

.75000E+C0 

•63760E-01 

.750005+00 

.809895-01 

•75000E+00 

• 11911 E+00 

•75000E+00 

. 7 50Q0E + 00 

. 19309E + 00 

. 75 OOOE + OO 

.236016+00 

1 75000E +00 

• 28 203E + 00 

.750005+00 

.330365+00 

•75C00E+00 

.75000E+00 

• A3077E+00 

. 75000E+00 

•A 81255+00 

.750005+00 

.530936+00 

• 75000E +00 

.57917E+00 

• 75000E+00 

•75000E + 00 

♦ 669A A E + 00 

• 75000E+00 

* 71090 E + 00 

.750005+00 

.7500CE+00 


J*1/K»1^X*1 

9 IMAX (X/Y/Z 

ARRAY) 







0. 

* 301 1 AE-01 

* 301 1 A E-01 

*68 190 E-0? 

• 3011 A E-01 

*3011 AE-01 

.136A0E-01 

. 3011AE-01 

•3011AE-01 

. 2QAQ0E-01 

. 3011A £-01 

. 30 11AE-01 

♦ 27370 E-01 

•3011A6-01 

•3011A6-01 

.3AA60E-01 

• 3011 AE-01 

•3011AE-01 

.A2120E-01 

♦ 301 1 A E-01 

• 3011 AE-01 

• 5 1190E-01 

. 3011 A£-01 

.3011AE-01 

• 63230E-01 

• 3011 AE-01 

.3011AE-01 

.80520E-01 

. 301 1 A E-01 

♦ 3011A E-01 

. 10580 E+00 

• 3011 AE-01 

• 30 11 AE-01 

.1A120E+00 

♦3011 fc E-01 

.3011AE-01 

. 17660E+00 

. 3011 AE-01 

. 3011AE-01 

.212006+00 

. 301 1 AE-01 

•3011 AE-01 

. 2 A7A0E + 00 

. 3011 AE-01 

. 3011AE-01 

. 2828 OE + 00 

. 3011 AE-01 

♦ 301 1 AE-01 

•31620E+00 

. 3011 AE-01 

*30 11 AE-01 

.353605+00 

. 3011 A6-01 

•3011AE-01 

.389005+00 

. 301 1 A E-01 

. 301 1 AE-01 

♦A2AAOE+00 

. 301 1 A£-01 

♦ 3011 A E-01 

• A59806 +00 

• 3011 AE-01 

.30 11 AE-01 

. A9520E +00 

•3011AE-01 

•3011*5-01 

•530605*00 

#3011*5-01 

•30114E-01 

• 56600E tOO 

•3011*E-01 

• 301HE-01 

.601A0E+00 

« 3011 AE-01 

• 3011AE-01 

.636006+00 

• 301 1 A E-01 

♦ 3011A E-01 

• 67220E +00 

• 3011 A B-01 

• 30UAE-01 

•70760E+00 

. 301 1 A E-01 

. 301 1 A£-01 

•7A300E+00 

. 301 1 AE-01 

.3011 AE-01 

•778A06+00 

. 3011 AE-01 

• 3011AE-01 


* .X73V2t+JJU 


•38021E+00 

•625ABE+00 


.277755-01 
•15393E+00 
*38021 E+00 
.625A7E+00 


•277755-01 

.1539AE+00 

.380215+00 

•625A7E+00 


•277755-01 
•1539AE+00 
• 3 6021 E+00 
•625A7E+00 


ITER NO.j IY 10CJ/ IYL0CK; YERRttX; 1ZLQC J> IZ LOCK; ZERRMX * 100 


15 10 


1669 A AE-OA 10 15 


16698 1E-0A 



*** CONVERGED SOLUTION ♦*+ 

CORNER FLOW VITH FIUET 
COORDINATE SYSTEM* S IZE-30X23X23 

JMAX • 23 KMAX • 23 

MAX ITERATIONS PERFORMED • 1X9 

MAX Y-ERROR • .998592E-05 A7 LOC. J>K « 15> 10 

MAX Z-ERRQR « .999905E-05 AT LOC# J>K ■ 10, 15 


J * 1 

.301146-01 
.41333E-02 
0# 

0# 

0# 

( Y-Z ARRAY) 

♦ 30114E-01 

• 73335 E— 01 
•218426+00 
f 4 44 05 E +00 
« 67204 E + 00 

• 26723E-01 
• A3390E-03 

0. 

0* 

0. 

•33673E-01 

#933796-01 

♦ 25933 E + 00 

• A 9200E +00 
•71209E+00 

.22417E-01 

0# 

0# 

0. 

0. 

J • 2 

•33673E-01 
« S 3075E-02 

• 3 562 IE-0 2 

• 3 560 46-0 Z 

• 35 603E-02 

( Y-Z ARRAY) 

# 267 23 E-01 

* 7 2769 E-01 
•218006+00 
• AA339E+00 

• 67201 E + 00 

•3170AE-01 

•A270CE-02 

• 356146-02 
•3060AE-02 

• 35603E-02 

• 3170 A E-01 
•92804 E-01 
•25397E+00 
.491876+00 
•71 200E+00 

• 27768E-01 
•360A7E-02 
• 35609E-02 

• 35603 E-02 
• 35603E-02 

J - 3 

t 36732E-0 1 
.1447GE-01 
•36290E-G2 

* 6 6 2 18E-0 2 

• 0 5 2 12 E-02 

(Y-Z ARRAY) 
.224176-01 

• 71835 E-01 

• 217 AO E + 00 
•443656+00 
•6719 7 E+00 

« 37 626E-01 
•93AA7E-C2 
•86257E-02 
• 66 215E-02 
•86212E-02 

*27768 E-01 
.919306-01 
•258A6E+00 
.A9169E+00 
•71206E+00 

•3A772E-01 

.87795E-02 

• 06239E-0 2 
.862146-02 

• 8621 2E-02 

J » A 

.46487E-01 

• 24 677E -0 1 
•16A21E-01 
•163G2E-01 

• 1 6 3 6 6c-Q 1 

(Y-Z ARRAY) 
•16C05E-01 
• 69891E-01 
• 2164 UE + OO 
• AA330E+00 
•67191E+00 

• 45759E-01 
.1390AE-01 

• 16 A06 6-0 1 
•16387E-01 

• 16 3 866-0 1 

• 2 10 75 E-01 

• 9 0370E-01 

• 2 5769 E + 00 

• 491 A 1E + 00 

• 7 1 2 03 E + 00 

• A3977E-01 

• 16950E-0 1 

• 1 6397 E-01 
« 16306 E-01 
•163C6E-01 

J « 5 

.57S20E-01 

• 3988 6E-Q 1 
#273916-01 
•277636-01 
• 27753E-0 1 

(Y-Z ARRAY) 

•10370E-01 

.658956-01 
• 21 51 1 E+00 
• AA279E+00 
•67182E+C0 

♦572160-01 

• 329 OQE-Ol 
•27631E-01 
•27759E-01 
» 277 53E-01 

•1MJ8E-01 

• 87 A9 2 E-01 

• 25655 E + 00 
.A9101E+00 
•71199E+00 

♦56Q48E*01 

• 29503E-01 
•27799E-01 
• 277 56E-01 
i 27753E-01 

J * 6 

• 733355-0 1 
•594996-01 

(Y-Z ARRAY) 
• 41333E-0* 
•59496E-01 

• 72 7695-01 

• 52 4 58E-01 

#830756-02 
• 8 27 0 A E-0 1 

• 71 835E-01 
• A750AE-01 


•38732E-01 
•11807E+00 
• 303036+00 
#539306+00 
•75000E+00 

♦16805E-01 

0# 

0# 

0# 

.46487E-01 
•14734E+00 
* 3488 6E + 00 
•58538E+00 

•10370E-01 

0. 

0. 

0# 

•57820E-01 
• 10092 E + 00 
•39610E+00 
•62974E+00 

•37626E-01 
•1 174CE+00 
♦30273E+00 
# 53920E + 00 
•75C00E+00 

•21875E-01 
•356646-02 
• 35607E-02 
♦35603E-02 

• 45759E-01 
•1467GE+00 
• 34861E+00 
• 58531E + 00 

• 14938 E-01 
•356386-02 • 
• 35605E-02 
•35603E-02 

• 57216 E-01 
• 18043E+00 

• 39590 E+00 
•629696+00 

• 3A772E-01 

• 11662 E + 00 

• 3023CE + 00 

• 53906E + 00 
*75000 E + 00 

• 29284E-01 
• 06566E-02 
•86228E-02 
•86213E-02 

* 43977E-01 

• 14599.E + 00 
.348266+00 
. 58521E+00 

• 21 831E-01 
#863616-02 
« 86 222E-02 
• 86212E-02 

• 56048E-01 
# 17973E+00 
• 39561E+00 

• 62962 E + 00 

• 29283E-01 
•11523E+00 
♦30166E+00 
•53B85E+00 
.750006+00 

• 39775E-01 
• 165A7E-01 
* 16393E-01 
♦ 163 86E-01 

* 39775E-01 

• 14474E + 00 
.347726+00 
• 58505E+00 

• 32741E-01 
.164546-01 
,163906-01 
♦16386E-01 

• 532 27E-01 

• 17864 E+00 
•395176+00 

• 62952 E+00 

•21831E-01 

•11296E+00 
. 30071E + 00 
#538556+00 
•75000E+00 

♦53Z27E-01 

• 2 8369E-01 
.27781E-01 
♦ 27755E-01 

♦327A1E-01 

•14203E+OO 

.34694E+00 

.584836+00 

♦47fc34E-01 

•28021E-01 
•27770E-01 
• 27754E-01 

• A 7 633 E-01 

• 17703 E + 00 

• 394 53 E+00 
•62937E+00 

•1447EE-01 
• 109A2E + Q0 

•69892E-01 

•45139E-01 

i2A67tE-01 

#140056+00. 

♦65895E-01 

• 39 8 86 E— 01 

. # 1 7 47I»C xf\t\ 



•A3793E-01 

• 21322E + 00 

• A3605E-01 

•25A98E+00 

• A3506E-01 

• A3395E-01 

•AA210E+00 

* A3382E-01 

• A90A7E + 00 

• A337 AE-01 

• ^ 3 3 6 5 £ — *0 1 

• 67169 E +00 

♦A3365E-01 

*71193 E + 00 

♦ A3365E-01 

J • 7 

• 93379E-0 1 

(Y-Z ARRAY) 
•A3390E-03 

♦ 9280 AE-01 

* A2700E-02 

• 9 1930 E-01 

•8270AE-01 

* 5 2A 57E-01 

•76576E-01 

• 765 75E-01 

•710A3E-01 

• 6A 6 6 8 E — G 1 

• 2 1077E + 00 

* 6A230E-01 

• 25296 E + 00 

• 63982 E-01 

• 63 702 E-01 

* 44122 E+00 

• 63669 E-01 

• A8979 E + 00 

•636A9E-01 

•63625E-01 

• 671 5 AE + Q0 

• 6362AE-01 

*711 86E+00 

•6362AE-01 

J » 8 

* 1 1 807E + 00 

(Y-Z ARRAY) 

0. 

•117A8E+00 

•360A7E-02 

•11662E+00 

• 109A2 E + 00 

• A750A E-01 

* 10 A58 E +00 

• 710 A2 E-01 

•99A21E-01 

• 90 663 E-01 

•20783E+00 

• 0999 AE-01 

• 2 5053E + 00 

.89A88E-01 

• 83897E-01 

. AA017E + 00 

• 888 27E-01 

♦A8896E+00 

•88783E-01 

•88732E-01 

•67136E+00 

• 88729E-01 

•71178E+00 

•88729E-01 

J * 9 

• 1A73AE +00 

C Y-2 ARRAY) 

0. 

• 1A676E+00 

•3568AE-02 

•1A599E+00 

♦1A005E+00 

•A5138E-01 

• 13629E+00 

• 67A60E-01 

•13193E+00 

» 12231E +00 

« 20A62E+00 

* 1209 2 E +00 

•2A783E+00 

•12006E+00 

• 1 1900E + 00 

* A 389 6 E + 00 

• 11 6 87E + C0 

• A 8802 E + 00 

• 11879E+00 

• 11 869E + 00 

♦67115E+00 

•11069E+00 

• 7 1 1 66 E + 00 

• 11869E+00 

J ■ 10 

• 18092E+00 

(Y-Z ARRAY) 
0* 

♦ 18 0 A3 E+00 

•35630E-02 

• 17973 E + 00 

« 17A75E +00 

• A A189E-01 

♦ 17175E + 00 

• 65595E-01 

• 16616E+00 

• 15853c + 00 

•20152E+00 

•15671E+00 

•2A508E+00 

•155A9E+00 

•153BAE+G0 

• A3766E+00 

•15363E+00 

• A8701 E + 00 

•153A9E+00 

• 1533 3 E +00 

•67092E+00 

« 15332E+00 

•71158E+00 

• 1 533 2E + 00 

J - 11 

t 218426 +00 

(Y-Z ARRAY) 
0* 

•21G00E+00 

•35621E-02 

« 217A0E + 00 

• 21323E +00 

• A3792E-01 

• 21070 E+00 

•6A607E-01 

•20783E+00 

f 198bfcE+00 

• 19686E+00 

• 19679E + 00 

* 2A2 55 E + 00 

•19528E+00 

•19306E+00 

* A3633 E + 00 

♦ 19275 E +00 

•A8596E+00 

• 1925 5 E + 00 

•19230E+00 

•67068E+00 

•1922QE+00 

♦711 A7E+00 

•19228E+00 

J * 12 

.259335+00 

(Y-Z ARRAY) 

Oi 

•25S97E+00 

*3B6H§-6§ 

•25gAfc§40G 

• 25A98E+00 

• A3605E-01 

•25296E+00 

•6A229E-01 

*25053 E+00 

• 2A256E + 00 

• 19679 E +00 

•2A0AAE+00 

* 2 A0A3E + 00 

•23B78E+00 

f 23610E+00 

• A3505E+0Q 

• 23570E + 00 

• A8A9AE + 00 

•235A3E+00 

1 23 5 10E + 00 

•670A5E+00 

* 23507E + 00 

*711 37E+00 

♦23507E+00 

0 • u 

•30303E+00 

(Y-Z ARRAY) 
0. • 

« 30272 E +00 

• 35609 E-02 

•30230E+00 

•299A1E+00 

.A3506E-01 

* 29775E +00 

•63981E-01. 

•29575E+00 

* 


1 



' • 


• 

* a * •* f vtr^vvr • - • 

•299A1E+00 
•5381 AE+OO 
•75000E+00 

• A3A50E-01 

• A3370E-01 

* *345888+00 
• 58A53E+00 

• 43 A16E-01 
• • A3367E-01 t 

•39367E+00 
.. .62917E + 00 

•98447E-02 

.msae+oo 

•29775E+00 

•53762E+00 

.75000E+00 

• 90371E-01 

* 67A61E-01 
•63839E-01 

• 63637E-01 

• 1890 AE— 01 
•13629E+00 
•3AA52E+00 
• 58A16E+00 

• 87A93E-01 
* 65596E-01 
• 6375AE-01 
•63629E-01 

•32980E-01 

• 17175 E+00 

• 392 57E + 00 
•62 891 E+00 

•87795E-02 
•99A19E-01 
• 29575E+00 
• 53700E+00 
•75000E+00 

• 11523E+00 
•95219E-01 
• 89190E-01 

• 8 8757E-01 

* 16950E-01 

• 1 3193E + 00 
• 3A288E+00 

• 58371 E + 00 

• 11296E + 00 

* 92 A72E-01 
• 89009E-01 

. • 887A1E-01 

.29503E-01 
• 16-816E+00 
.39124E+00 
•62861E+00 

•86566E-02 

.95218E-01 

•29350E+00 

•53629E+00 

.75000E+00 

•1AA7AE+00 * 
• 1277QE+00 
• 1195AE + 00 
•1187AE+00 

♦ 165 A7E-01 
• 12777E+00 

* 3A102E + 00 

* 5 8 3 1 9 E ♦ 0 0 

♦14283E+00 
•12452E+00 
♦ 11921E+00 
•11871E+00 

•28369E-01 

• 16443 E+00 

• 3 8973E+00 
•62826E+00 

.86361E-02 

•92471E-01 

•29114E+00 

♦53552E+00 

♦75000E+00 

•1786AE+00 
• 16 A A3E + 00 
« 1 5A69E+00 
•153A1 E+00 

•16A5AE-01 
•12A52E+00 
♦3390AE+00 
• 5826 A E+00 

* 17703E+00 
• 1611 1E+00 
•15A18E+00 
• 15336E+00 

• 28021E-01 
• 161UE+00 
*38811 E+00 
•62768E+00 

♦86290E-02 

.90882E-01 

•28887E+00 

.53A72F+00 

.75000E+00 

•216A8E+00 
•20A62E+00 
•19A2AE+00 
• 192A2E+00 

•16A21E-01 
•1 2231E+00 
• 33709E+00 
• 58206E + 00 

* 21 511E + 00 
•20152E+00 
• 19353E+00 
• 1923AE+00 

• 27891E-01 
• 1 58 53E + 00 
• 30647E+00 
•627A9E+00 

i86257§“02 

•89993E-01 
.28686E+00 
• 53399E+00 
•75000E+00 

•25769E+00 

• 2A783E+00 
• 23756E + 00 
•235Z6E+00 

•164065*01 

•12092E+00 
• 33528E+00 
* 581 A9E + 00 

•256555+00 

•24508E+00 
* • 23670E+00 
•23516E+00 

•270J15-O1 

•15670E+00 

•36492E+00 

;;.62710E+00 

.86239E-02 

«8948EE-01 

• 30166E+00 
♦29350E+00 

• 1 6397E-01 
• 12006E + 00 ... 

• 30071E+00 
... *29HAF+ftn 

•27799E-01 



•28887E+00 

• 19526E+00 
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CORNER FLOW WITH FILLET 
COORDINATE SYSTEM.SIZE 30X 23X 23 

— INITIAL GUESS — 

JMAX 23 EM AX 23 
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CORNER FLOW WITH EILLET 
COORDINATE SYSTEM. SIZE 30X23X23 

— CONVERGED SOLUTION — 

JM AX 23 EH AX 23 



III. Programmer/Analyst Section 


The program is written such that the array size in n, C 
direction can be varied from problem to problem by the use of 
test editor on CDC 6600 machine. The array sizes are coded to 
be IMXLL, JMXLL, KMXLL in £ , n and ? direction respectively. 
These are changed to the appropriate numeric number by the test 
editor prior to compiling the program for execution. The job 
control cards, input explanation and other helpful hints are 
provided towards the beginning of main program COORDC. The 
program is extensively commented. Though the computer code 
is not structured in the true sense, extreme care was taken 
to avoid branching back in the code unless absolutely necessary. 
A brief discussion of each subroutine and their function follow. 


A. Main Program COORDC 

This program performs the following function, 

i) Initialize variables to 0.0. 

ii) Read in input data. If the field size is greater than the 
maximum field size allowed then the program aborts with a 
diagnostic message. The maximum field size allowed is 
ILIMIT, JLIMIT, KLIMIT and the field size input is IMAX, 
JMAX, KMAX . ILIMIT, JLIMIT, and KLIMIT are set to numerical 
value by the text editor. The storage requirement in N-S 
code require that IMAX, JMAX, KMAX be equal to ILIMIT, 
JLIMIT, KLIMIT respectively. 

iii) Specify initial guess. 

iv) Check to see if initial guess is to be printed, plotted or 
not. If so, then the appropriate subroutines are called. 

v) Compute attraction parameter Q and R. 
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vi) Start computing solution. Do loop for maximum iterations 
(ITERMX) is set. For each iteration subroutine CALCOR 

is called. CALCOR computes solution for each iteration 
and also passes back maximum error and its location. 

vii) Print maximum error and partially converged solution 
every 50 iterations. Partially converged solution is 
printed only if the flag to print the initial guess is 
set. 

viii) Check maximum error in y and z for convergence. If con- 
vergence is reached, then an exit from the DO loop is 
taken. 

ix) After ITERMX iterations or convergence, which ever comes 
first, some more print/plot and storing of data is done. 

x) Check to see if partially or converged solution is to be 
printed/plotted . 

xi) If convergence has occurred, then store solution on unit 1 
If convergence has not occurred, then the maximum errors 
in y and z computation at the last iteration and their 
location are printed out. 


B. Subroutine CALCOR 

This subroutine advances the solution 1 iteration. The 
maximum errors and their location are computed and passed back 
to the calling program COORDC through the argument list. The 
argument list is. 


YERRMX - 

maximum 

error in y 

ZERRMX - 

maximum 

error in z 

IYLOCJ - 

maximum 

y-error 

J location 

IYLOCK - 

maximum 

y-error 

K location 

IZLOCJ - 

maximum 

z-error 

J location 

I Z LOCK - 

maximum 

z-error 

K location 



c. 


Subroutine IPRTC 


This subroutine is used to print initial guess, partially 
converged or converged solution. The argument list is r 
ISOLN = 1 indicates initial guess 

= 2 indicates partially or fully converged solution 
KKK = applicable only if ISOLN =2 

= 0 indicates partially converged solution . 

= 1 indicates fully converged solution. 

YERRMX, ZERRMX, IYLOCJ, IYLOCK, IZLOCJ, IZLOCK = same as in 
CALCOR. 

D. Subroutine IPLTC 

This subroutine is used to plot initial guess, partially 
converged solution or converged solution. The argument list 
is same as IPRTC. 

E. Block Diagram 

A Block Diagram is presented on the next page. 
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